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Abstract 

We consider random Schrodinger equations on Z d for d > 3 with identically dis- 
tributed random potential. Denote by A the coupling constant and ipt the solution with 
initial data tpQ. The space and time variables scale 

A -2-«/2 )t ^ A -2-K with 

< k < Ko(d). We prove that, in the limit A — > 0, the expectation of the Wigner 
distribution of ipt converges weakly to a solution of a heat equation in the space variable 
x for arbitrary L 2 initial data. The diffusion coefficient is uniquely determined by the 
kinetic energy associated to the momentum v. 

This work is an extension to the lattice case of our previous result in the continuum 
[8], [9]. Due to the non-convexity of the level surfaces of the dispersion relation, the 
estimates of several Feynman graphs are more involved. 
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1 Introduction 



We consider the time evolution of the Anderson model [2 J given by the random Schrodinger 
equation 

id t Mx) = HMx), ^e£ 2 (Z d ), t€R (1.1) 
on the d-dimensional square lattice Z d . The Hamiltonian is given by 

H = -±A + XV U , (1.2) 
where A > is the coupling constant. The kinetic energy operator on £ 2 (7j d ) is given by 

(Af)(x):=2df(x)- £/(x + e) (1.3) 

|e|=l 

and the random potential is given by 

V u {x) = V M v *( x ) ■■= v * S x , a (1.4) 

a£Z d 

where {v a : a E Z d } are i.i.d. random variables and 6 x , a = 1{ X = a } is the usual 
Kronecker delta function and 1{-} is the characteristic function. We will work in d = 3 
dimension, but our results and proofs extend to any d > 3 in a straight-forward manner. 

We study the long time evolution of the equation (jl.ip . For a large coupling constant, 
A > Ao, the spectrum of H is almost surely pure point and the dynamics is localized 
[2], [TTJ [T| . It is conjectured, but not yet proven, that the spectrum is absolutely continuous 
and the dynamics is diffusive if A < Ao is sufficiently small. We will investigate the 
dynamics in this regime in the scaling limit, when time diverges as A —* 0. 

Up to time scales t ~ A -2 the dynamics is kinetic, typically with a finite number 
of collisions. The evolution on macroscopic space scales, x ~ A~ 2 , is given by a linear 
Boltzmann equation [131 El El [12] . As the long time limit of the Boltzmann equation 
is the heat equation, the quantum evolution on scales t S> A -2 is expected to exhibit 
diffusive behavior. In this paper we prove this statement up to time scale t ~ A~ 2 ~ K 
with a positive k. 

We have proved the same statement (with a somewhat bigger k) for a random 
Schrodinger operator in the continuum, R rf [U [9]. The history of this problem and 
related works are summarized in [8] and will not be repeated here. 

Our reasons to extend this work to the lattice case are: (i) to show that the methods 
initiated in [6l[5], and later developed in [HI [9] for longer time scales, can be applied to 
a lattice setting as well; (ii) to make a connection with the extended state conjecture of 
the Anderson model. 
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Anderson (de)localization is a large distance phenomenon, thus no physical difference 
is expected between the lattice and continuum models. The localization proofs, however, 
are typically simpler in the lattice models because of certain technical difficulties due to 
the ultraviolet regime in the continuum model. 

We now explain briefly the differences between the continuum and lattice models in 
our analysis of the derealization regime. The finite momentum space is an advantage in 
this regime as well; the artificial large momentum cutoff's introduced for the continuum 
model in [3] are not necessary here. Moreover, the computation of the main term is 
more direct since the Boltzmann collision kernel is homogeneous on the energy shells 
(compare (I2.19P below with (2.19) of |8j). In particular, the diffusion coefficient can be 
computed explicitly. 

However, an important technical estimate is considerably more involved for the lattice 
case. The complication stems from the non-convexity of the isoenergy surfaces of the 
lattice Laplacian. The isoenergy surfaces are the level sets, S a := {p G T rf : e(p) = a}, 
of the dispersion relation 



Our approach heavily uses estimates that integrals of resolvent functions, \a— e(p) + 
W7| -1 , T) <C 1 are concentrated near E Q . In particular, our most involved estimate (Four 
Denominator Lemma [3.4p translates naturally to a specific decay property of the Fourier 
transform of measures supported on the hypersurfaces S Q . Such decay bounds are readily 
available for surfaces with non-vanishing Gauss curvature, but we were unable to find in 
the literature the necessary estimate for surfaces with Gauss curvature vanishing along 
a submanifold. Since in the energy range a £ (2, 2d — 2), the Gaussian curvature of 
the level sets T, a vanish along a codimension one submanifold (see Figure [1]), we had to 
prove this estimate separately [10] . Another related bound (Two Denominator Lemma 
17. 7p will be proven in this paper. 

The analogous bounds in the continuum model, with dispersion relation e(p) = \p 2 , 
are much easier. The actual proofs (e.g. Proposition 2.3 in [4J or Lemma A.l in [8]) use 
the explicit form of e(p), but the key fact is that the level sets of the dispersion relation 
are convex. 

Most importantly, the Two-Denominator Lemma is much stronger in the continuum 
case; the estimate in Lemma 17.71 carries a diverging factor ^~ 3 / 4 ~ K . The corresponding 
bound (formula (10.28) in [8j) is only logarithmically divergent in d = 3. This makes 
the estimates of the error terms in [9] much easier. 

We mention that the difficulty related to the non-convexity of the level sets is present 
in the proofs of the Boltzmann equation from a lattice model (3j [12] as well. On the 
kinetic time scale, however, weaker bounds (called crossing estimates) were sufficient 
(Lemma 7.8 of [3] or Assumption (DR4) of p2]). 
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Figure 1: Level set of e(p) = a for 2 < a < 4 
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The main inputs we use from our other works are: (i) the stopping rules (Section 3 
of [2]); (h) the basic integration procedure for non-repetition Feynman graphs (Sections 
8-10 of [8]); (iii) the organization of the estimates on the error terms (Section 4 of 
[9]); and (iv) the Four Denominator Lemma proved in |10j . The new ingredients are: 
(i) the estimates of the error terms (with the Four Denominator Lemma replacing the 
strong form of the Two Denominator Lemma not available in the lattice case); (ii) the 
computation of the main term and the explicit diffusion coefficient. 

2 Statement of the main result 
2.1 Notations 

In this paper we consider the random Schrodinger operator (|1.2p acting on £ 2 (Z d ). The 
kinetic energy operator is given by the discrete Laplacian (|1.3p and V w (x) is the random 
potential (jl.4p . Denote the moments of the single-site random potential by =EcJ. 
We assume that mi = 771,3 = m 5 = 0, 771-6 < 00 > an d 7712 = 1 by normalization. Universal 
constants will be denoted by C and their value may vary from line to line. 
On the lattice (5Z) d , 5 > 0, we introduce the notation 



On the dual space (T/5) , with T := [—5, 5]) the integration refers to the usual Lebesgue 
integral 



We will use these formulas mostly for 5 = 1 in which case we will drop the subscripts 
(5'L) d and (T/5) d indicating the integration domains in (|2.ip . (|2.2p . The letters x,y,z 
will always denote lattice variables, while p, q, r, u, v, w will be reserved for d-dimensional 
momentum variables on the torus. This notation will distinguish between the two inte- 
grations (|2.ip and (|2.2p . Note that as 5 —* 0, both integrals converge to the standard 
Lebesgue integral on M. d . 

For any / £ l 2 ((5'L) d ) the Fourier transform is given by 




(2.1) 




(2.2) 




(2.3) 



where p = (p^ 



. 6 (T/5) d . The inverse Fourier transform is given by 
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The notations and (•) will be used according to convenience. For functions denned 
on the phase space, f(x,v), with x G v G (T/5) d , the Fourier transform will 

always be taken only in the space variable, i.e. 

M, v) := / e- 2m ^f(x, v)dx, £ G (T/5) d . 

We also remark that addition and subtraction of momenta will always be defined on the 
torus, i.e. with periodic boundary conditions. 

The Fourier transform of the kinetic energy operator (jl.3|) on ^ 2 (Z d ) is given by 

(FAf)(p) = -2e(p)/(p) , 

where e(p) is the dispersion law defined in ([13]) . For /i : T d C and an energy value 
e G [0, 2cZ] we introduce the notation 

[h)(e) := J h(v)S(e - e(v))dv : = jf h(q) (2.4) 

where dz^(g) = dv e (q) is the restriction of the d-dimensional Lebesgue measure to the 
level surface 

S e := {q : e{q) = e} C T d . 
By the co-area formula it holds that 

i-2d 



[h](e)de = / h(v)dv . (2.5) 
o J 

We define the projection to the energy space of the free Laplacian by 

( h(v) ) e := ^4 - wh ere $(e) := [11(e) = fs(e- e(u))du . (2.6) 
*(e) ' J 

Define the Wigner transform of a function -0 G £ 2 (1j d ) 

W4x,v) := 2 d £ e 2 ™^Wmz), * G G T rf . (2.7) 



2 

T/ + 2 — 2a; 

Notice that the position variable x takes values on the refined lattice (Z,/2) d . It is easy 
to compute the Fourier transform of W^(x,v) in x using (12. 3D with 5 = 1/2, and one 
obtains 

%(M =$(y- |)^(« + |) , (e(2T)SeT d . (2.8) 
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Note that W^(£, v) is periodic with respect to the double torus (2T) d . It is easy to check 
that W^(x,v) reproduces the correct marginals in the following sense 

u, ( \a J A^)\ 2 if x g z d 



r 

J(Z/2) d 



W4x,v)dx = 2- d Wtp(x,v) = $(v)\ 2 , (2.9) 



ase(2/2) 



and in particular 

/ / Wjf,(x, v)dvdx - 

J(Z/2) d J 

Recall that momentum integrations on unspecified domains are always considered on 
T d . Define the rescaled Wigner distribution as 

W$(X, V) := e~ d W^ (—, , X £ (eZ/2) d , V G T d . (2.10) 
Its Fourier transform in X is given by 



W^,V)=T{W^-,V))(C) = 4>(v- £ 4)^(v+ £ 4) , S, G (2T/e) d , V G T d 



To test the rescaled Wigner transform against a Schwarz function 0(x,v) on IR rf x T rf , 
we introduce the scalar products: 



{0,W$):= I dX [ dvO(X,v)W^(X,v) 

J(eZ/2) d J 

( 3, W$ : = [ d£ / dv 6(S, v)W°(eS, v) . (2.11) 

J(2T/e) d J 

By unitarity of the Fourier transform, ( O, Wl) = { O, W^). 

One may, of course, also define the Wigner transform of a function tp G £ 2 (1i d ) by 
first defining W^(£, v) for £ G (2T) d by ([22]), then taking the inverse Fourier transform 
in £ to reproduce the lattice function defined in (12. 7p . An interpretation of the Wigner 
function as a distribution starting from (|2.8p for any £ G M rf was given in |12j . 

2.2 Main theorem 

The weak coupling limit is defined by the following scaling: 

T:=et, X:=ex, e = A 2 . (2.12) 
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It was proved [3] that in the limit e — > the Wigner distribution W^ T {X, V) converges 
weakly to the solution of the Boltzmann equation 

c\F T {X, V) + sin(27ry) • V X F T (X, V) 

= 2tt / dU5(e(U) - e(V)) F r (X, U) - F T (X, V) , (2.13) 

with velocity vector sin(2-7ry) := (sin(27rFi), sin(27rV^), sin(27rV3)) = ^pVe(V). Here 
Ft(X, V) is the time dependent limiting phase space density with X £R d ,V £ T d . Note 
that the Boltzmann equation can be viewed as the generator of a Markovian semigroup 
on phase space. In particular, all correlation effects become negligible in this scaling 
limit. 

In this paper we consider the long time scaling 

x = \- K / 2 ~ 2 X = e~ l X, t = \- K ~ 2 T = s- 1 \- K/2 T, s = \ k/2+2 (2.14) 
with k > 0. Our main result is the following theorem. 

Theorem 2.1 Let d = 3 and ifto G £ 2 (lj d ) be an initial wave function. Let ip(t) = 
solve the Schrddinger equation hl.l\) . Let 0(x,v) be a Schwarz function on R d x T d . 
For almost all energies e € [0, 2d], [IV'o^)! 2 ] ( e ) is finite and for these energies let f be 
the solution to the heat equation 

8 T f(T, X, e) = V x • D(e)V x f(T, X, e) (2.15) 

with the initial condition 

f(0,X,e) := 5(X)[\Mv)\ 2 ](e) 

and the diffusion matrix D 

( sin(27ruW)-sin(27n;0'))) 
DM:= K ^ i,j = 1,2,3. (2.16) 

Then for k < hq and e and A related by (|2.14p . the Wigner distribution satisfies 

lim / dX fdv 0(X, v)BW" (x . K . 2T) (X, v) = [ dX [dv 0(X, v)f(T, X, e(v)) , 

(2-17) 

and the limit is uniform on T £ [0, To] with any fixed Tq. In d = 3 dimensions, one can 
choose kq = 1/9800. 
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Remarks, (i) We stated the Theorem and will carry out the proof in d = 3 for 
simplicity but our method works for any d > 3. 

(ii) The analogous theorem for the continuum model was proved in [8] and [9] with 
a somewhat larger threshold for k. The threshold is obtained from technical estimates, 
it has no physical relevance and it can be improved with a more careful analysis. 

(iii) By the symmetry of the measure (-) e under each sign flip Vj —* — vj and by the 
permutational symmetry of the coordinate axes, we see that D(e) is a constant times 
the identity matrix: 

(sin 2 (2W 1 ))) 
D ij (e)=D e 6 ij , D e := X 



2?r $(e) 

in particular we see that the diffusion is nondegenerate. 

(iv) The diffusion matrix can also be obtained from the long time limit of the Boltz- 
mann equation f)2. 13j) (see also the explanation to Figure 1 in [8]). For any fixed energy 
e, let 



L, 



J{v) := J du a(u,v)[f(u) - /(«)], e(v) = e , (2.18) 

be the generator of the momentum jump process on S e with the uniform stationary 
measure (-) e . The collision kernel is 

a(U, V) := 2vr5(e(C/) - e(V)). (2.19) 

The diffusion matrix in general is given by the velocity autocorrelation function 

Dij(e) = \ dt (frin(27r«W(*)) • sin(27ro^(0))) , (2.20) 

where v(t) is the process generated by L e . Since for any fixed V the collision kernel 
cr(U, V), is symmetric on the energy surface in the U variable, the correlation between 
v(t) and v(0) vanishes after the first jump and we obtain (|2,16j) . by using 



J dua(u,v) = 27r<3?(e) , e(v) 



e . 



3 Preparations 
3.1 Renormalization 

The purpose of this procedure is to include immediate recollisions with the same obstacle 
into the propagator itself. Without renormalization, these graphs are exponentially large 
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("divergent"), but their sum is finite. Renormalization removes this instability and the 
analysis of the resulting Feynman graphs will become simpler. 

The self-energy operator is given by the multiplication operator in momentum space 

9{p) := e(e(p)), 9(a) := lim 9 e (a) , (3.1) 

where 

9.(a) := / %— . (3.2) 

J a - e(q) + is 

The existence of the limit and related properties of are proved in Lemma lA.ll 
We rewrite the Hamiltonian as 

H = H + V, 

where 

H :=u{p) := e(p) + \ 2 9{p), V := XV - \ 2 9{p) . (3.3) 

We note that our renormalization is only an approximation to the standard self-consistent 
renormalization given by the solution to the equation 

w (p) = e (p) + A 2 lim / . . dg . . . (3.4) 

V ' ' J L0{p) - u(q) +IE K ' 

Due to our truncation procedure, the definition (|3.1|) is sufficient and is more convenient 
for us. 

We also need a few properties of the function u(p). First, uj(p) is symmetric under 
the permutation of the coordinate axes. Second, e(p) is symmetric for reflections onto 
any coordinate axis: e(pW, . . . = e(±p^\ . . . , ±pW) for any choices of the signs 

and these properties are inherited by 9(p) and w(p). Moreover, e(| — p) = 2d — e(p), 

1 TnoroTnrci iHllnyl — — C-^iOrl — rv \ anH 

2 • V2' 2' " " " » 1> 



where \ := (i, |, . . . , i). Therefore 0(a) = — 0(2d — a) and 



1 



^2 -?) = -%)■ (3-5) 
These relations allow us to restrict our attention to the subdomain 

V:={p=(p( 1 \p ( - 2 \...,pW)£T d : 0<pV<pM<...<pW ,^> (i) <|} (3.6) 

i 

of the momentum space. On this domain we have the following estimate on 9(p): 
Lemma 3.1 For any d>3 there exist universal positive constants c\,oi such that 

\9(p)\ < c 2 , (3.7) 

- c 2 \p\ d - 2 < Im 9(p) < - Cl \p\ d - 2 (3.8) 

for any p £ D. 
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Using this bound and the symmetry properties of ui{p) above, we easily arrive at 
Corollary 3.2 For d > 3 we have 

Imu(p) < -c 3 X 2 [D(p)] d ~ 2 , p G T d (3.9) 

for some C3 > 0, where 

D(p) : = min ||p — v \ : v = or v = ^ rfc — , ± — , . . . , } 

is i/ie distance between p and the set consisting of the origin and the vertices ofT d . □ 

We remark that this estimate would produce logarithmic corrections in d = 2 dimensions, 
this is one of the reasons why the proof is somewhat simpler in d > 3. 

Proof of Lemma \3.1\, By the co-area formula (|2.5p and recalling the definition (|2.6p 
we can write e (a) as 

QJa)= r $(s), with *(«)= / -^44 • (3-10) 

V 7 io a-s + K 7e s |Ve(g)| V 7 

Because of the symmetries of e(q), it is sufficient to study 

du(q) 



$ si := 



s a n© |Ve(g)| ' 



i.e. where the integration is restricted to T>. 

The critical values of e(q) are the even integers, 2m, between and d, and within V 
they correspond to the critical points 

( 1 1 

p m = 0,0, . . 




m times 

with m < ci/2. If s is away from a neighborhood of the critical values of e(q), then <&(s) 
is a bounded function with a strictly positive lower bound. 

If s is near a critical value, s = 2m + /3 with /3 sufficiently small, then we can write 

where Xm is a smooth cutoff function around the critical point p m . The second integral 
is a regular function in s that is separated away from 0. The first integral can be brought 
into the following normal form by a smooth local coordinate transformation: 



Jc(3) \Q\ 



C(f3) 
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where % is a smooth cutoff function around and 

d—m d 

C{(5) :={ 9 £l d : ti U) f ~ E [« (J) ] 2 = P } ■ 

j=l j=d—m+l 

It is a straightforward calculation to see the following behavior of the function > 
for small (5 

$#(/?) = O^i" 1 ), ^#((3) = 0^^ if m = (3.11) 
= O(l), —$#(/?) = O(l) if l<m<d/2. 

From these estimates, $ is differentiable away from 0, bounded everywhere and 
behaves as s 1 / 2 near in (i = 3 dimensions. For higher dimensions is bounded. From 
the formula (|3.10p . we can rewrite as 

f a ds 

e e (a)=/ <J?( a - S ). (3.12) 

Using the above properties of $ we can take the limit e — > and define 0(a) := 
lim e ^o+o &e(a). If we write 0(e) = 1Z(e) — iI(e), where 1Z(e) andX(e) are real functions, 
and recall Im(x + 10) _1 = — ir5(x), we have 

Im0(a) = -7T J 5(e(q) - a)dq = -vr$(a) . (3.13) 

Using the properties of we can also check that for any < e < d 

$( e ) = /(e)ef- 1 , K{e) = a(e)(l + el -1 ), J(e) = 6(e)e^ 1 (3.14) 

where a,b,f are bounded functions, uniformly separated away from zero and /(e) > 0. 
We also have 1(a) = 7r<&(a). The estimates in Lemma [3.11 then follow from the fact that 
e(p) > cp 2 on T>. Later we will also need the bounds 

|*'(e)| <C(l + ei- 2 ) |0'(e)| < \K' (e)\ + \1' (e)\ <C(l + ef~ 2 ) , (3.15) 

that can be proven by a similar analysis. □ 

The following lemma collects some estimates on the renormalized propagators we 
shall use to prove Theorem 12.11 Its proof will be given in the Appendix. These are 
technical bounds and their meanings will become clear when they are used. 
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Lemma 3.3 Suppose that X 2 > n > \ 2 + 4k with k < 1/12. Then we have, 

su p / 1 TT-r-i ^ c \ lo § A l ( 3 - 16 ) 

a J \a - u(jp) + irj\ 



<C a \- 2 ^- a K (3.17) 



and for < a < 1 

f dp 
sup / | — — - 

a J \a - U(p) + I7]\ z a 

For a = 0, the following more precise estimate holds. There exists a universal constant 
Cq such that 



■I 



sup I , ^.f . 2 ■ i • r„A' 

\a — uJ(p) — irj\ 



[ ! N A ' d f 7 , r < l + CoA 1 - 12 - 

■J \ot 



(3.18) 



a,p,r J I" - w(p + r) - «7| |/? - w(p - r) + ir/| 

3.2 Truncation 

For any real number a we define 

|||a||| := min{|a|, \a — 2|, |a — 3|, \a — 4|, |a — 6|} (3.19) 

in the (i = 3 dimensional model. The values 0, 2,4, 6 are the critical values of e(p). The 
value a = 3 is special, for which the level surface {e(p) = 3} has a flat point. For |||a||| 
separated away from zero, the following key estimate holds (the proof is given in [10J): 

Lemma 3.4 [Four Denominator Lemma] For any A > rj there exists C\ such that for 
any a G [0,6] with |||a||| > A, 

dpdqdr ^ ^ M _ 14 



a — e(p) + irj\ \a — e(q) + irj\ \a — e(r) + irj\ \a — e(p — q + r + u) + irj\ 
uniformly in u. □ 



< C A | log 77I 



In general, in d > 3 dimensions, |||a||| is the minimum of |a — d\ and of all \a — 2m\, 
0<m<d. 

We will prove the main Theorem 12.11 under the assumption that the initial wave 
function in Fourier space is continuously differentiable, HV'ollc 1 < 00 j an( i ^ satisfies the 
following condition: There exists A > such that 

> 3A on the support of Vo • (3.20) 
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Once the theorem is proven for such initial data, we can easily extend it for the 
general case. Since e(p) is a Morse function, for any positive A > we can define a 
smooth cutoff function < x A < 1 on T d with the property that 

|||e(p)||| > 3A on the support of x A 

and 

hin y |1 - x A (p)\ 2 dp = . 

Moreover, we assume that x A ^ X A pointwise if A < A'. The wave function will be 
decomposed as 

where ipi , ip2 are defined in Fourier space as 

Mt) := e- UH ( X A ^o) Mt) := e~ UH [(l - X A )4>o] ■ 

The Wigner transform enjoys the following continuity property: if the random wave 
function is decomposed as ifi = ifi\ + ip2 > then 

E{d,W%)-E(6,WD < ( [ sup|a(^)|d£yE[||V>i|| 2 + ||V 2 || 2 ] -Ell^P 

(3.21) 

(see Section 2.1 of [6], but due to a misprint, the HV^H 2 term was erroneously omitted). 
Since 

II^WH = \\e~ im [(l - X k )U\\2 = ||(1 - X A )^o|| 2 , 
by monotone convergence we see that 



lim 

A^O 



E(S,W$ w )-E<a,W£ 



(t)> 



uniformly in t (and thus in e). This means that the truncation procedure is continuous 
on the left hand side of ([2TH) . 

Similarly, on the side of the limiting heat equation, we can define f A (T,X, e) to be the 
solution to (I2TT5D with initial data / A (0, X, e) := 5(X) [|x A ^o| 2 ] (e). Clearly [|x A ^ | 2 ] (e) 
monotonically converges to [|V>o| 2 ] (e) as A — > for any e such that |||e||| ^ 0. Therefore 
/ A (0, X, e) converges to /(0, X, e) in L^{dX), and thus the same statement holds for the 
solution 

f A (T,X,e)^ f(T,X,e) in L\dX) 

for almost all e and uniformly in T. But then the right hand side of (|2.17p is also 
continuous as A — ► 0. 

The condition, HV'ollc 1 < 00 can 8 ^ so be removed by an analogous truncation argu- 
ment, see Section 3.2 of [8] for details. 
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4 The proof of the Main Theorem 



2.1 



The structure of the proof is the same in the continuous and in the lattice case, so this 
section is almost identical to Sections 4 and 5 of [§]• There are three minor differences 
in the structure. First, in [8] the single site potentials were indexed by the finite set 
{1,2,..., M}, while here they are indexed by Z d (|1.4p . Second, in the continuous case 
the problem was first restricted to a finite box (Section 3.3 [8]) and this restriction 
was removed at the end of the analysis. This complication is absent here. Finally, the 
Boltzmann collision kernel contains an additional factor in the continuous case (compare 
(|2.19p with the corresponding formula (2.19) of [8]). 

We expand the unitary kernel of H = Hq + V (see (I3.3P ) by the Duhamel formula. 
For any fixed integer N > 1 

N-l 

i/H := e~ itH ^ = Mt) + ^N(t) , (4.1) 

n=0 

with ^ 

ip n (t) := (-i) n f [dsj]^ 1 e- iSn+lHo Ve- iSnHo V ...Ve- islHo ^ (4.2) 
Jo 

being the fully expanded terms and 

:= H) f ds e-^ H V^ N ^{s) (4.3) 
J o 

is the non-fully expanded or error term. We used the shorthand notation 



ft ft ft n n 

/ «:= / •••/ (n d ^>(*-E- 

Jo Jo Jo K jJi ' v jr x 



Since each potential V in (|4.2p . (j4.3j) is a summation itself, V = —X 2 6(p) + J2 a 
both of these terms in (|4,2p and (|4.3p are actually big summations over so-called elemen- 
tary wavefunctions, which are characterized by their collision history, i.e. by a sequence 
of obstacles and, occasionally, by 0(p). Denote by T n , n < oo, the set of sequences 

7 = (7l,72,--.,7n), j 3 £Z d U{$} (4.4) 
and by the associated potential 



Wzf if 76 Z rf 

■\ 2 9{p) if 7 = 1?. 



The tilde refers to the fact that the additional {■!?} symbol is also allowed. An element 
7 G Z d U {1?} is identified with the potential and it is called potential label if 7 G 
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otherwise it is called a i^-label. Potential labels carry a factor A, -^-labels carry a factor 
A 2 . 

For any 7 6 T n we define the following fully expanded wavefunction with truncation 

VH7 := H)"" 1 A da ;]i W^e-^W^ . . . e-^W^e-^ 110 ^ (4.5) 
Jo 

and without truncation 
^ : = ['[dsj]^ 1 e- i *+ lHo WT n e- i ' nHo W^_ 1 . . . e-^^W^e-** 1 ^^ • (4.6) 

JO 

In the notation the star (*) will always refer to truncated functions. Note that 

V>t, 7 = H) [ dse'^-^As^ ■ 
Jo 

Every term (|4.6p along the expansion procedure is characterized by its order n and 
by a sequence 7 £ T n . The main term is given by the sequences that contain different 
potential labels only. Their set is defined as 

Tf :={ 7 = ( 7l ,..., 7fe ) : lj € Z d , 7i ^ 7j - if i} • (4-7) 

Let 

denote the corresponding elementary wave functions. 

The typical number of collisions up to time t is of order X 2 t. To allow us for some 
room, we set 

K := [X~ 5 (X 2 t)} , (4.8) 

([ • ] denotes integer part), where 5 = 5(k) > is a small positive number to be fixed 
later on. K will serve as an upper threshold for the number of collisions in the expansion. 

The proof of the Main Theorem 12.11 is divided into three theorems. The first one 
states that all terms other than < k < K, are negligible: 

Theorem 4.1 (L 2 -estimate of the error terms) Let t = 0(A~ 2-K ) and K given by 
(|4.8p . If K < Ko(d) and 5 is sufficiently small (depending only on k), then 

K-i 2 

hmE||v^-E^|| =°- 

~* k=0 

In d = 3 dimensions, one can choose kq(3) = qskq. 
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The second theorem gives an explicit formula for the main terms, ip?l., expressed in 
terms of the ladder diagrams. We introduce the renormalized propagator 

1 

R v {a,v) : 



a — u(v) + ir] 



Theorem 4.2 (Only the ladder diagram contributes) Let k < e = A 2+K / 2 , 
rj = X 2+K , t = 0(A~ 2_ft ), and K given by (|4.8j) . For a sufficiently small positive 5 and 
for any 1 < k < K we have 



nmw 2 



V x {t,k) + 0[\^ 9K -°^ (4.9) 
(0 5 EW|»r) = W x {t,k,0) + o(\ts- 9k -° {5 A (4.10) 



as A -C 1. Here 

Vx(fk) :=-^r/ J dad/3 e^-^y (jl^i) l^(Pi)| 2 
fc+i 

w A (t,A;,o) :=^^ y j_ oo dadi3ei{a ~ l3)t j {2T/e)d d ^ j (n dv a(e ' ufe+i) ^° (e ' ui) 

fe+l 7— 

^^("'"i + fW^-f)- ( 4 - 12 ) 

i=i 

We adopt the notation 0(<5) in the exponent of A. This always means (const.) 5 with 
universal, explicitly computable positive constants that depend on n and that can be 
easily computed from the proof. We note that the definition (|4.12p does not apply 
literally to the free evolution term k = 0; this term is defined separately: 

W x (t, k = 0, 0) := f d£ f dv e it£ ^ e{v) e 2txHva6 ^ 6{£, v)W^ (eC, v) . (4.13) 

J{2T/e) d J 

The third theorem identifies the limit of W\(t, k, O) as A — > with the solution 
to the heat equation. 
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Theorem 4.3 (The ladder diagram converges to the heat equation) Under the 
conditions of Theorem \4-S\ and t = \~ 2 ~ K T, we have 

K-l 

Urn W x (t, k, O) = / dX / dv 0(X, v)f(T, X, e(v)) , (4.14) 
fc=o ^ ^ 

where f is the solution to the heat equation (|2.15p . 

Proof of the Main Theorem \2.1\ using Theorems \4-l\ \4-%\ an d \4-3\ We compute the 
expectation of the rescaled Wigner transform, ~EW£ = EWI, tested against a Schwarz 
function, (O, EWf) (see (|2.1ip ). By combining Theorem 14.11 with the L 2 -continuity 
of the Wigner transform (|3.2ip . it is sufficient to compute the Wigner transform of 
tp(t,K) := Ylk=o $tk- ^ ne Wigner transform Wm^K) contains a double summation 

K-l 

w i>(t,K)= E ^(••- )«(•••)• 

fe,fc'=0 

The potential labels are not repeated within ip and if). Moreover, the expectation of a 
single potential in (|4.6p is zero. Thus the potential labels in the ip and ip must pair, in 
particular taking expectation reduces this double sum to a single sum over k 

K-l 
k=0 

By using (gJDJ) and (|4TT4]I together with K = 0(\~ K ~ S ), we obtain Theorem^ □ 



5 Stopping rules 

The Duhamel expansion (|4.ip allows for the flexibility at every expansion step N to 
decide if the full evolution e -*(*- s )# j n (|4,3p is expanded further or not. The decision is 
based upon the collision history of the expanded terms. The stopping rules organize the 
expansion. The basic idea is to expand up to the identification of the main terms, but 
not to expand error terms unnecessarily further. The stopping rules are identical in the 
continuous and lattice cases and they were given in Section 3 [9] in full details. Here we 
only summarize the concepts informally and refer to [9] for the precise definitions. 

In a sequence 7 E T n we identify the immediate recollisions inductively starting from 
71 (due to their graphical picture, they are also called gates). The gates must involve 
potential labels and not Any potential label which does not belong to a gate will 
be called skeleton label. The index j of a skeleton label 7j in 7 is called skeleton index. 
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1 2 3 4 ' ! 5 6 7 ' ! 8 9 10 11 1 ! 12 13 

O Non-skeleton index — Q — 0— Gate 

i 

1 i 

# Skeleton index 



Figure 2: Skeleton indices and gates for 7 = (a, a, 6, 6, c, d, d, #, e, e, /). 



The set of skeleton indices is S^j) . The terms are never part of the skeleton. This 
definition is recursive so we can identify skeleton indices successively along the expansion 
procedure. Notice that the last skeleton index may become a gate index at the next step. 
Fig. [2] shows an example for these concepts, the formal definition is given in Definition 
3.1 [9j. The number of skeleton indices will be denoted by ^(7) := |5(7)|. Let £(7) 
denote the number of ^-labels in 7. Recalling that potential labels carry a factor A and 
^-labels carry a factor of A 2 , we let 

r(l) ■= ^[n-k(nf)]+t(fr) (5.1) 

denote the total A 2 -power collected from non-skeleton indices. 

Sequences where the only repetitions in potential labels occur within the gates are 
called non-repetitive sequences. A special case of them are the sequences in T^ r (j4.7|) 
that contain no gates or ^-labels. The repetitive sequences are divided into the following 
categories. If two non-neighboring skeleton labels coincide, then the collision history 
includes a genuine (non-immediate) recollision. If a skeleton label coincides with a gate 
label, then we have a triple collision of the same obstacle. If two neighboring skeleton 
labels coincide and are not immediate recollisions because there are gates or 0's in 
between, then we have a nest. The precise definitions are given in Definition 3.3 [9]. 
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Genuine recollision Triple collision with a gate 

a b b a 

-Q — © — •- 



# Skeleton index 

O Non-skeleton index 



Nest 



We stop the expansion at an elementary truncated wavefunction (|4.5p characterized 
by 7, if any of the following happens: 

• The number of skeleton indices in 7 reaches K. We denote the sum of the truncated 
elementary non-repetitive wave functions up to time s with at most one A 2 power from 
the non-skeleton indices or 0's and with K skeleton indices by £ ■ The superscript 
(< 1) refers the number of collected A 2 powers. 

• We have collected A 4 from non-skeleton labels. We denoted by tp^Jf ast the sum 
of the truncated elementary wave functions up to time s with two A 2 power from the 
non-skeleton indices (the word last indicates that the last A power was collected at the 
last collision). 

• We observe a repeated skeleton label, i.e., a recollision or a nest. The corresponding 
wave functions are denoted by ^' reC > V^f 1 f }' nest . 

• We observe three identical potential labels, i.e., a triple collision. The corresponding 
wave functions are denoted by Tp[~ 1 *• 

Finally, denotes the sum of non-repetitive elementary wavefunctions without 

truncation (i.e. up to time t) with at most one A 2 power from the non-skeleton indices or 
9's and with k skeleton indices. In particular, the non-repetition wave functions without 
gates or •& (denoted by ip^k aD ove) contribute to this sum. For convenience, we will 

rename ip < fl' nr '■= V^fc indicating explicitly the number of A 2 -powers from non-skeleton 
labels, r = 0. 

This stopping rule gives rise to the following representation that is intuitively clear. 
The detailed proof is given in Proposition 3.2 of [9]. 
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Proposition 5.1 [Duhamel formula] For any K > 1 we have 

^ = e- itH ^ = ^4f' nr (5.2) 

fc=0 

• Ca -i(t-s)H ( i(<i),nr . ST^ ( , {2), last , , (<l),rec . ,(<l),nest, , (<l),*ri\ 1 
~ t J dSe Y*s,K + Z^{^*s'k +^*~s,k + ^*l,k +Ks,k )}■ 

6 Graphical representation 

For Theorems 14.11 and 14.21 we need to compute expectations of quadratic functionals 
of elementary wavefunctions t/Jt,j- This computation is organized by Feynman graphs. 
In the continuum model, the computation was shown for the non-repetition terms in 
details in Section 6 of [8] and the precise definitions of the Feynman graphs were given 
in Section 7 of [SJ. For the lattice model, the Feynman graphs are very similar but 
somewhat simpler because the ultraviolet regime is absent and the Boltzmann collision 
kernel (|2.19p depends only on the energy. We will not repeat the complete arguments 
here, but we introduce the necessary modifications for the lattice case. 

6.1 Circle graphs and their values 

We start with an oriented circle graph with two distinguished vertices denoted by 0, 0*. 
The number of vertices is N. The vertex set is V, the set of oriented edges is C{V). For 
v G V we use the notation v — 1 and v + 1 for the vertex right before and after v in 
the circular ordering. We also denote e v _ = (v — l,v) and e v+ = (v,v + 1) the edge 
right before and after the vertex v, respectively. For each e G C(Cn) we introduce a 
momentum w e and a real number a e associated to this edge. The collection of momenta 
is denoted by w = {w e : e G £(V)} and dw = ® e dw e is the Lebesgue measure. The 
notation v ~ e indicates that the edge e and adjacent to the vertex v. 
Let P = {P fl : H G 1} be a partition of the set V \ {0, 0*} 

v\{o,o*} = Up m , 

(all are pairwise disjoint and non-empty), where I = I(P) is the index set to label the 
nonempty sets in the partition. Let m(P) := |/(P)|. We will call the sets P-lumps 
or simply lumps. Two elements of v, v' G V are called P-equivalent if v, v' belong to the 
same lump of P. 
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We will assign a variable (called auxiliary momenta), £ T d , jjL £ I(P), to each 
lump. We will always assume that the auxiliary momenta add up to 

M6/(P) 

The set of all partitions of the vertex set V \ {0, 0*} is denoted by Py. For any P C V 
we let 

L + (P) := {(v,v + 1) £ £(V) : « + 1 P, u £ P} 

denote the set of edges that go out of P, with respect to the orientation of the circle 
graph, and similarly L_(P) denote the set of edges that go into P. We set L(P) := 
L + (P)UL_(P). 

For any £ £ M. d , define the following product of delta functions 



A € (P,w,u) :=5^+ Y ±w -) II 6 { E ±™e-^J, ( 6 - 2 ) 

eeL±({0*}) M6/(P) eeL±(P M ) 

where u := {u^ : // £ /(P)} £ T rf is a set of auxiliary momenta. The sign ± indicates 
that momenta w e is added or subtracted depending whether the edge e is outgoing or 
incoming. We also recall that all momentum variables live on the torus, in particular 
addition is also defined on T d . The dependence on £ will be mostly omitted from the 
notation; A = A^. All estimates will be uniform in £. 

Summing up all arguments of these delta functions and using (16. If) we see that these 
delta functions force the two momenta corresponding to the two edges adjacent to to 
differ by £: w e — w e > = £ for e £ L + ({0}), e' £ L_({0}). This will correspond to the 
momentum shift in the Wigner transform (|2.8p . 

With these notations, we define for any P £ Py and r/ > the P-value of the 
partition 

E(P,u,cx):=\ N - 2 [ J] 1 ■ | A(P,w,u)|^oKo + )ll4Ko-)l • (6.3) 

The prefactor A^ -2 is due to the fact that in the applications all but two distinguished 
vertices, 0, 0*, will carry a factor A. The E- value depends on the parameters A and rj but 
this is omitted from the notation. In the applications, the regularization will be mostly 
chosen as rj = A 2+K . 

We will also need a slight modification of this definition, indicated by a lower star in 
the notation: 



P*(P, u, a) := \ N ~ 2 f I] dw e T] , 1 A(P,w,u)|^oKo + )ll^oK JI • 

(6.4) 



ee£(V) eG£(V) 
^£({0*}) 
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The only difference is that the denominators carrying the momenta associated to edges 
that are adjacent to 0* are not present in E* v . We will call the truncation of E. We 
will see that Feynman diagrams arising from the perturbation expansion can naturally 
be estimated by quantities of the form (|6.3p or (|6.4p . 

The formulas (|6.3p and (|6.4p are the lattice analogues of (7.5) and (7.6) of [8] in the 
continuous model, but the momentum cutoffs, the polynomially decaying factors for the 
special set Q and the non-trivial collision kernel B are absent. 

Following the continuum model, we also define four operations on a partition given 
on the vertex set of a circle graph and we estimate how the the E- value changes. The 
estimates are somewhat simpler for the lattice case, so here we just summarize the results 
and prove an additional estimate (see (|6.8p below). For the details see Lemma 9.5 of [8] 
and Appendix C of [9]. 

Lemma 6.1 (Operation I. Breaking up lumps) Given P = {P^ : jjL G /(P)} G 

T\, we break up one of the lumps into two smaller nonempty lumps; P v = P u i UP„" with 
P v iV\P v ii = §. Let P* = {P u i , P v u , P^ : /i£/(P)\{i/}} denote the new partition. Then 

£ w (P,u,a)< J drE^(P*,u*(r,u),a) , (6.5) 

where the new set of momenta u* = u*(r, u) is given by it* := u^, fi G L(P) \ {v} and 
u*, = u v — r, u*„ = r. In particular, 

sup £7 W (P, u, a) < sup E (if) (P* , u, a) . □ (6.6) 

u u 

The notation E^ simultaneously refers to E and i.e. to formulas with and without 
truncation. 

Lemma 6.2 (Operation II. Removal of a single vertex) Let v E V \ {0,0*} be a 
vertex and let P G Vv such that P a = {v} for some a £ L(P), i.e. the single element set 
{v} is a lump. DefineV* := V\{v}, £(V*) := C(V)U{(v-l,v+l)}\{(v-l,v), (u,v+l)}, 
i.e. we simply remove the vertex v from the graph and connect the vertices v — l,v + 1. 
Let P* := P \ { {v} } G V v * , then 

supE(^(P,u, a) < Xr]^ 1 sup £'(„) (P* ,u* ,a) . (6.7) 

u u* 

Furthermore, if v and 0* are neighbors in the graph, then we have the following 
stronger estimate for the truncated value: 

sup-E*(P,u,a) < A (3d+ |a e |) sup£(P*, u*, a) , (6.8) 

u u* 
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where e is the edge connecting 0* with its neighbor other than v. 

Finally, if both neighbors of 0* , v ^ v' , form single lumps in P, then both of these 
lumps can be simultaneously removed to obtain a partition P* := P \ {{w},{f'}} with 
the estimate 

supK»(P,u,a) < A 2 sup£(P*,u*,a) . (6.9) 

u u* 

Proof. Estimates (16. 7ft and (16. 9ft were proven in Appendix C of [9j. For the proof of 
(|6.8p . notice that if v + 1 or v — 1 is 0*, then the denominator \a e — uj(w e ) +irj\~ 1 is not 
present in the truncated value for e = (v — 1, v) or e = (v, v + 1). Therefore the proof 
of (|6.7p can be repeated without paying the price. The extra factor estimates the 
missing denominator \a e — oj{w e ) + irj\ < \a e \ + 3d. □ 

Lemma 6.3 (Operation III. Removal of half of a gate) Let v, v + 1 G V \ {0, 0*} 

form a gate in a partition P G Vy, i.e. P a = {v,v + 1} for some a G I(P)- Define 
V* :=V\{v + 1}, £(V*) := £(V) U {(v,v + 2)} \ + l),(u + l,v + 2)}, i.e. we 

simply remove the vertex v + 1 from the circle graph with the adjacent edges and add a 
new edge between the vertices v,v + 2. Let P* G "Py* ^ e P after simply replacing the 
lump {v,v + 1} with {v}. Then 

£7 w (P,u,a) < A|log77| S w (P*,u,a) . □ 

Lemma 6.4 (Operation IV.: Removal of a gate) Let v, v + 1 G V \ {0, 0*} /orm a 
gaie m P G Vy, i.e. P a = {v,v + 1} for some a G I(P). Define V* := V \ + 1}, 
£(V*) := J C(V)U{(u-l,« + 2)}\{(u-l,t;),(u J i; + l),(u + l,« + 2)}, i.e. we simply 
remove the gate. Let P* G Py* ^ e P after removing the lump {v,v + 1}. Then 

sup-E(*)(P,u, a) < \ 2 r]~ 1 \ logry| sup£ , (jf )(P*, u*, a) . □ 

u u* 



6.2 Feynman graphs 

Feynman graphs are special circle graphs that naturally arise in the perturbation expan- 
sion. Consider the cyclically ordered set V„ )Tl ' := {0, 1, 2, . . . , n, 0*, h', n' — 1, . . . , 1} and 
view this as the vertex set of an oriented circle graph on N = n + n' + 2 vertices. We 
set L n := {1, 2, . . . n} and I n > := {1, 2, . . . , n'} and the vertex set can be identified with 
V n ,n' = In U J n ' U {0, 0*}. The set of edges C(V n>n >) is partitioned into C(V n>n ') = C U C 
such that £ contains the edges between I n U {0,0*} and C contains the edges between 

/n'U{0,0*}. 
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Let V n n > be the set of partitions P of the set I n U I n >. Let G = G(P) be the set 
of edges that enter a single lump and let g(P) := |G(P)| be its cardinality. In case of 
n = n', we will use the shorter notation V n = V nj „, V n = V n ,n etc. In our applications 
we always have |n — n'\ < 4 and n, n' < K. 

Let Q : M. d — > C be an arbitrary function that will represent the momentum depen- 
dence of the observable. For convenience, we can assume that ||Q||oo < 1- For a,/?€l, 
P G T ) n,n'- > and we define 

tt(P, a, /?) := \ n + n '+a( p ) fdwTT -L _ TT _L _ (6.10) 

xA^(P,W,U = 0) Yl [-6( w e)} II {- e ( W e)]^( w e 0+ )^(We -)Q ^(We ,_+W eot+ ) , 

eecnG eeCDG 

with w := {w e : e G C U £}. The subscript £ will mostly be omitted. 

The truncated version, V^(P, a, (3), is defined analogously but those a and /? denom- 
inators are removed that correspond to e G L({0*}). We set 

2tr? /• /-4<Z 

V W (P) := — ^ y y dad/3^ a -^V w (P,a,l3) (6.11) 

and 

' " a) 



S w (P,u) := ^2 y y ^ dad/3 S w (P,u, 



where a in £V*) (P, u, a) is defined as a e = a for e G £ and a e := /3 for e G £. We 
will call these numbers the F-value and E- value of the partition P, or sometimes, of the 
corresponding Feynman graph. Strictly speaking, they depend on the vector £ and on the 
function Q as well; when necessary, we will use the notations V^^P, Q), E^(P, u, Q), 
etc. 

Clearly 

|V W (P)| < (CA)f( p ) £7 w (P,u = 0) . (6.12) 

If u = 0, we will use the notation Eu\(P) := EuU~P, u = 0). 

As we will see, in the graphical representation of the Duhamel expansion what we 
really need is 

V&(P) := ^0 |^dad/?e^^V w (P,a,/3) , (6.13) 

i.e. a version of VL)(P) with unrestricted dad/3 integration (the circle superscript will 
refer to the unrestricted integration). However, the difference is negligible even after 
summing them up for all partitions. 
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Figure 3: Vertex and edge labels 



Lemma 6.5 Assuming that rj > A 2+4k and 1 < n + n' < 2K , we have 

Y, |V W (P) - V5j(P)| = 0(Al^')) . (6.14) 

T/ie same result holds if V^^P) were defined by restricting the a, [5-integral to any 
domain that contains [— 4d, 4d] x [— 4d, 4d]. 

Proof. Outside of the regime |a|, < 4d, at least either the denominators with a 
or with (3 in (|6.10p are uniformly bounded since < 2d + 1 for small A. The other 

denominators can be integrated out at the expense of (C|logA|) max ( n ' n ') +1 by using 
(|3.16p . The contribution to V^{P) from the complement of |ct|,|/3| < 4d is therefore 
bounded by 



i/ w (p)-y ( : } (p) 



< yi+n'+g(P)fQt log ^|\max(n,n')+l 



if A is small. Since the total number of partitions, \V n ^ n '\i is bounded by (n + n') n+n ' 
and n + n' <2K = 0(A" K ~ 5 ), we obtain (IBTT4"]) . □ 

Sometimes we will use a numerical labelling of the edges, see Fig, O In this case, we 
label the edge between (j — 1, j) by ej, the edge between (j, j — 1) by e~-. At the special 
vertices 0,0* we denote the edges as follows: e n+ \ := (n, 0*), e^TTJ '■= (0*,n'), e\ = (0, 1) 

and ej := (1,0). Therefore the edge set C = C(V n ^ n i) is identified with the index set 
I n+ i U I n '+i and we set pj := w ej , := w e ~. These two notations will sometimes be 
used in parallel. Note that we always have 

Pi-Pi=£- (6-15) 



6.3 Non-repetition Feynman graphs 

Let A n be the set of partitions of I n := {1,2,..., n}, i.e. A = {A^ : fi G -f(A)} G ,4 n 
if U^ e /(A)^4^ = -^n and the elements of A are disjoint. The sets in the partition are 
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labelled by the index set 1(A) and let m(A) = |/(A)| denote the number of elements in 
A. The elements of the partition A will also be called lumps. A lump is trivial if it has 
only one element. The trivial partition, where every lump is trivial, is denoted by Aq. 

A partition P G V n of I n U I n is called even if for any P^ G P we have \P„ H I n \ = 
\Pfi H I n \. In particular, in an even partition there are no single lumps, G(P) = 0. 

Let & n be the set of permutations on I n and let id be the identity permutation. Note 
that A G A n and a G © n , uniquely determine an even partition in P(A, a) G V n , by 
/(P) := 1(A) and P M := U a(A^). Conversely, given an even partition P G V n , we 
can define its projection onto I n , A := vr(P) G A n , by 1(A) := I(P) and A^ := P a n J n . 
We let 

6„.(P):={(re6„ : P(tt(P), a) = P} 

be the set of permutations that are compatible with a given even partition P. In other 
words, a G S n (P) if for each i 6 I n the pair (i,a(i)) belongs to the same P-lump. 
Clearly 

ie„(p)i= n (^)!= n ( 6 - 16 ) 

M e/(P) Me/(7r(A)) 

We will use the notation 

y w (A,a,Q) := y w (P(A,a),Q) (6.17) 

and similarly for E^ g and In the proofs, Q will be omitted. We also introduce 

c(A):= [J C (KI). ( 6 - 18 ) 
^eJ(A) 

where c(n) are the coefficients of the connected graph formula defined in (|6.23[) . With 
these notations we can state the representation of the non-repetition terms as a summa- 
tion over Feynman diagrams. 

Proposition 6.6 With Q = 1 and £ = we have 

mZkf= E E c(A)V£o(A,a,Q = l) (6.19) 

anc? with Q^(v) := ioe /ioue 

B(d,W; r )=J2 E < A ) / de^(A,a,Q c ). (6.20) 

aG6 fc A6^ fc ^(2T/ £ )^ 
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The proof is essentially given in Section 6 and Proposition 7.2 of [8] with a few minor 
modifications. The restriction to the finite box Al is absent. The expectation value of 
the potentials in the expansion of E||?/j"£|| 2 is given by 

k k 

(6.21) 

instead of (6.4) of [8]. The following identity was proven in Lemma 6.1 [8]: 
Lemma 6.7 For any fixed k, 

k 

c(n):= £ (-l)l r l, (6.23) 

rcA'„ 

r connected 

where K n denotes the complete graph on n vertices and \T\ denotes the number of edges 
in the subgraph T. The following estimate holds 

\c(n)\ < n n ~ 2 . □ (6.24) 



7 Proof of Theorem 4.2 



We recall that the E- and ^-values of the partitions depend on the parameters A, t, £ 
and k; a fact that is not explicitly included in the notation. In Sections [7] and [8] we will 
always assume the following relations 



v = A 2+K , t = \- 2 - R T, T E [0, T ], K = [A" 5 (A^)], k < K, (7.1) 

with a sufficiently small positive 5 > that is independent of A but depends on k. All 
estimates will be uniform in £ and in T E [0, To]. 



7.1 Estimates on graphs with high degree 

We recall the key definition of the degree of a permutation a E from Definition 8.3 
of [8]. Let a act on = {1, 2, . . . , k} and let a be its extension to {0, 1, . . . , k + 1}, by 
a(0) = 0, a{k + 1) = k + 1, otherwise a{i) = a{i). An index j E / := {1, 2, ... k} is called 
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Figure 4: Ladder graph 

ladder index of a if a(j) — 1 € {a(j — 1), a(j + 1)}. Let It be the set of ladder indices 
and i = £(a) := Finally, the degree of a is defined as 

deg(cr) :=k — £(a) . 

Starting with (|6.19p , we notice that the graph with the trivial partition Ao and with 
the identity permutation on 1^ gives the main term in Theorem 14.21 since 

V\(t,k) = V^ =0 (Aq, id) , c(Ao) = l. 

This graph is called the ladder graph (Fig. H|) . 

To prove that all other graphs are negligible, we first replace V°(- ■ ■ ) with V(- • • ); the 
error is negligible by Lemma [6 .51 We then first estimate V( A, a) for the trivial partition 
A = Ao, where every lump has one element. Since |F(A,cj)| < E(A, a, u = 0), the 
following bound is the key estimate: 

Theorem 7.1 Assume (|7.1|) with k < ^ and let a £ 6^. Then the E-value of the 
graph of the trivial partition with permutation a is estimated by 

sup# w (Ao,<7,u) < c(A^- 4K ) degM |lo-gA| 2 (7.2) 

u ' ' 

i/A< 1. 

The proof is essentially given in Section 10 of [5] with a few modifications that we will ex- 
plain in Section [7T21 below. This theorem is complemented by the following combinatorial 
lemma that was proved in [8] (Lemma 8.5): 

Lemma 7.2 Let k < 0(\~ K ~ S ), D > integer, and let 7 > k + 5 be fixed. Then 

A 7d0g(CT) < C»(a d(7 - k - 5) ) (7.3) 

CTgSj, 

deg(tr)>D 

for A«l. □ 

Since deg(o") > 2 if a 7^ id, from Theorem l7.1l and Lemma [7.2l we immediately obtain: 
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Proposition 7.3 Assuming (|7.ip with k < ^ we have 

^2 \V(A ,a)\ < o(\*- WK -°^ (7.4) 

for A < 1. □ 

For the general case A ^ Ao, we need to recall the notion of joint degree of a 
permutation a E and a partition A E Aj~ from Definition 9.1 of [8]. Let S'(A) be 
the union of non-trivial lumps in A, and let s(A) := |5(A)| be its cardinality. The joint 
degree of a pair (a, A) is given by 

q(a,A) := maxjdeg(cr), ^s( A )} ■ 

The following statement was proved in Lemma 9.3 [8] (recall the definition of the com- 
patibility and that of the projection tt(P) from Section f6.3|) . 

Lemma 7.4 For any even partition P E Vk there exists a compatible permutation a = 
a(P) E e fc (P) such that 

deg(a) > ~ S (tt(P)) . □ (7.5) 

The following corollary shows that the estimate of a general partition can be reduced 
to that of a trivial partition with the help of Lemma EH The proof is somewhat simpler 
than Corollary 9.4 of [S] in the continuous case. 

Corollary 7.5 Given a E 6k and A E Ah, we have, for k < 

sup£ w (A,a,u) < C|logA| 2 (A^- 4K ) 9(A,,T) . (7.6) 

u ^ ' 

Proof of Corollary |7.5[ We define a permutation a* := a* (A, a) as a* := a if 
deg(cr) > |s(A), and a* := <t(P(A,ct)) otherwise. By Lemma IT~4l we have deg(c*) = 
q(A,a). Clearly P(A,<t) = P(A,cr*), in particular, (A, a, u) = £7*)(A, er*, u). By 
Operation I. we can artificially break up all non-trivial lumps in A and use the auxiliary 
momenta to account for the additional Kirchoff rules. Using the estimate ()6.5[) and that 
J dr = 1, we immediately see that 

sup (A, a* , u) < sup £(*) ( A , a*, u) , 

u u 

and Theorem 17.11 completes the proof. □. 

Finally, we have the following bound on the summation of general non-repetition 
graphs. The proof is the same as the proof of Proposition 9.2 of [H] but the estimate 
(|7.6p replaces Corollary 9.4 of [8] in the argument: 
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Proposition 7.6 We assume (|7.ip . Let D > and s > 2 be given integers, let q := 
max{L>, t^s}. For any k < jj^ and 5 < 5(k) sufficiently small, we have 

]T Yl sup^ W (A,«7,u)|c(A)| < c(A^- 9K -°W) 9 |l°gA| 2 . □ (7.7) 

deg(<j)>D s(A)>s 

Finally, the proof of Theorem 14.21 follows from Proposition 17.61 exactly as the proof 
of Theorem 5.2 in the continuous case explained in Section 7.3 of [8] (Proposition 16.61 
replaces Proposition 7.2 of |8j and Lemma 16.51 replaces Lemma 7.1 of [3]). □ 



7.2 Proof of Theorem O 

The proof of Theorem 17.11 follows the integration scheme presented in Section 10 of 
[8] with a few modifications. Apart from the finite momentum space, the finite cutoff 
in the a, (5 variables and the simpler collision kernel, the main difference is that the 
Two Denominator Lemma is weaker in the lattice case (compare Lemma [7. 71 below with 
Lemma 10.5 in [S]). This results in weaker A-exponents in the estimates and eventually 
a smaller threshold kq in the main result. 
We define 

||| ||| :=77 + min{| 9 - 7 0')|, j = 1, 2, . . . , 2 d } (7.8) 

with 7«) = la(j) and a(j) = (ai(j'), d2(j), ■ ■ ■ , o-d{j)) is the dyadic expansion of j — 1. 
In other words, |||g||| — r\ is the minimal distance of q from the critical points of e(p) 
(measured on the torus T d ). This is not a norm, but it satisfies the triangle inequality, 

lb + < Ibl + IIMI- 

For any index set I' C I := {1, 2, . . . , k + 1}, any |7'| x (k + 1) matrix M and any 
vector b = &2, . . . , bk+i) & we define 

E(I',M,b): = A 2fc sup [f dad/3 sup f TT dp,- ^ f 



X 



II ,. -mi .. „...., .... r-r l- < 7 - 9 > 



u(pi) - ir}\ \/3 -oo(Y,j=i M ijPj + "i) + H 

where u € M 3 is an additional dummy momentum and 

b • p := hipi + 6 2 P2 + • • • + b k+1 p k+1 G M 3 . 

We will use the notation E{I' , M, 0) defined exactly as (|7.9p but without the factor 
|||b • p + ^ HI 1 in the integrand and without the supremum over v. We will refer to this 
case as choosing the "empty vector" b = 0. Notice that the definition (|7.9p is somewhat 
simpler than the corresponding (10.14) of [8]. 
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For any permutation a G &k we associate a (k + 1) x (k + 1) matrix M(a) according 
to (8.7) of [8]. This matrix encodes the momentum dependencies in the delta function 
A in the E- value of the partition P(Aq, a). We have 



sup£(Ao,a,u) < \\^ \\i o E(I,M(a)^) . (7.10) 

u 

The truncated version, E 1 *, can be estimated by the untruncated one E*(. . .) < CE(. . .) 
because in the regime \a\, \(3\ < Ad every propagator is bounded from below. 

An easy estimate is available for Eu\(l , M(<r),0) by first separating all but one a 
and (3 denominator by a Schwarz inequality (k > 1): 

E(I,M(a),$) <A 2fc sup [dp [ 1^ , ^ (7.11) 

a J J J -id l« - v(pi) - vr)\\P - u(pi +ux) + 

k+1 ^ k+1 



(: ;2 l« - u{pj) - ii]\ 2 ' jJ- 2 \a - uj(qj) + irj\ 2 

where qi = J2j Mij( a )Pj + The estimate for E* has one less denominator. Since 
M(a) is invertible with determinant ±1 (Proposition 8.2 of [8]), the contributions of 
the two terms in the square bracket are identical. To estimate the first term, we can 
integrate out all pj variables, j = 2, 3, . . . , k + 1, by using (|3.18p . yielding a factor 
(1 + CoA 1 " 12 ^)^ = 0(1) since k < K < A~ 1+12k . After integrating a, (3 and finally pi, 
we obtain 

sup E(I, M(a), 0) <C| log A| 2 (7.12) 

crGSfe 

sup E m (I, M((r),0) <CA 2 |logA| 2 . 

CT6Sfe 

Note that the squared denominators can be integrated out in arbitrary order, unlike in 
the continuum case (Section 10.1.2 of [8]). 

To obtain a bound of the form \(c°nst.)deg(a) for E(I,M(a),$), one has to gain a 
A-power from the non-ladder variables. This requires a successive integration procedure 
described in Section 10.3 of [8]. We will not repeat here the formal procedure, but just 
mention the basic idea. The difficulty is that for a general a each variable pj appears 
in many denominators in (|7.9p . To break this complicated dependence structure, a set 
of carefully selected /3-denominators in (|7.9p are estimated trivially by rj . Performing 
then the pj-integrations in a specific order, the remaining (3 denominators can be inte- 
grated out (together with the a-denominators) without losing further ry _1 factors. Each 
integration involves only two propagators (Two Denominator Lemma IT.Tj) . 
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In principle, if the propagators corresponding to the ladder indices are estimated by 
a Schwarz inequality argument (|7.1ip . one should gain A 2 from each non-ladder index. 
This would give a bound of order \ 2de s( a ) i n (|7.2p . modulo logarithmic corrections. 
Unfortunately, point singularities may arise from the repeated application of Lemma 
17.71 this necessitates the factor |||b • p + wfl -1 in (|7.9p . To avoid the accumulation of 
point singularities during the integration procedure, we estimate trivially not only the 
/3-denominators with ladder indices, but several other ones as well. This accounts for 
the smaller power in (|7.2p . 

The integration procedure removes the /^-denominators one by one. Each /^-denominator 
in (17. 9p is labelled by an index i G I and we will treat them in increasing order of the 
index i. The index set / = {1, 2, . . . , k + 1} is partitioned into six disjoint subsets, 

/ = I p U I v U h U I cs U I uc U I last (7.13) 

described in Definition 8.3 and Definition 10.3 of [8j. 

To bookkeep the integrations, in Section 10.3 of [8] we defined a sequence of matrices, 
a sequence of index sets, I^ h \ and a sequence of vectors, f»w, for h = 1, 2, . . . , k+1. 
We set E(h) := E(I^ h \ M^, b^ 1 )). Initially = /, AfW = M(» and t>(°) = 0, so 
from (THUD 

sup^(A ,a,u) < H^ollL^(O) • (7.14) 

u 

As h = 1,2,... increases, in each step we estimate E(h — 1) in terms of E(h). The 
estimate depends on the set where h £ I falls into according to the partition (|7.13p . 
The actual estimates are somewhat different in the lattice case than the corresponding 
bounds (10.19), (10.20), (10.25), (10.29), (10.32) and (10.35) of [8J in the continuous 
case. We will list the results only, the proofs are analogous to the arguments in [8]. 
In Case 1, h G I p , we have 

E(h - 1) < r]~ l E(h) . (7.15) 

In Case 2, h G Ie, we consider h, h + 1, . . . ,h + t — 1 € I? a maximal sequence of 
consecutive ladder indices (i.e. h — 1, h + r Ii with some r > 1), then 

E(h - 1) < C\- 2r E(h + r - 1) . (7.16) 

The proof is easier here in the lattice case: since the form factor B is absent in (j7.9[) . 
each ladder index can be integrated out independently by using (|3.18p . Note that the 
constant in (I7.16P is independent of r. 

In Case 3, h G I us , we use ()3.16|) and (|A.29P from the Appendix, i.e., the lattice 
versions of (10.23) and (10.24) from [8]: 

E(h - 1) < CV W2 | logrifE(h) < Cr ] - 1 - K E(h) . (7.17) 

In Case 4, h G I cs , we need the following lattice version of Lemma 10.5 in [8] that 
will be proved in the Appendix: 
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Lemma 7.7 [Two Denominator Lemma] For r\ = A 2+K we have 

1 1 , Cri- 7 /*~ K 

iup / ! =r\ — 

|a|,|/?|<4d 

Without the point singularity we have 



/l 1 Cv~ 7 ^ 8 ~ K 
i =T~\ r TTR 7 T \ i • i i if d P ^ O ( 7 - 18 ) 
|q - Lo(jp) - irj\ |/3 - uj(p + q) + ir]\ \jp - r\j \jq\l 



sup r dp < g^ 3/4 ' K (7 19) 

M,[/9[<4ei./ l« - w(p) - ir?| |/3 - a/(p + g) + i»7| _ |||g||| ' '"" 



Finally we have 



" da f __ » . _ < ^V^I . (7 .20) 



I a — u>(p) — ir/\ \a — co(p + q) + ir/| 



Notice that these bounds are weaker than the ones given in Lemma 10.5 [8] which 
themselves are not optimal. For example, the factor ^ _7 /8- K [ n (17.18P can be improved 
to ?7 -1 / 2 in the analogous estimate for the continuum model. 

Using this bound and following the argument of Case 4 in Section 10.3 of |8j, we have 

E(h - 1) < Ci]~ 7/8 ~ K E(h) . (7.21) 

In Case 5, h £ I v , we use 

f f dpdp' 1 

|a|,|/3|<4d r P J J \a - uj(p) - irj\ \a - To(p') - irj\ |/3 - u(p -p' + q) + ir]\ \\p - r\\ 

< Cr]- 7 / 8 - 2 * 

instead of (10.31) of [8]. This inequality follows from (|7.18p and (|A.29p . The same bound 
holds if the point singularity is of the form \\p ±p' — 7"|| _1 (compare with (10.32) of [S]) 
or if there is no point singularity at all. Following the argument of Case 5 in [8], we 
obtain 

E{h - 1) < Crf 7t8 - 2K E{h) . (7.22) 

if h e I v . 

In the last step, h = k + 1 E hast, we can estimate E(k + 1) by directly integrating 
out a and (3 similarly to (10.35) of [3] to obtain 

E(k + 1) < CA 2fc | log/yl 2 . (7.23) 

Combining fTTIjh ([7TTK|) - (f717]l and (17T2T1 - (17331) . and using that Case 2 has been 
applied not more than k — £ = deg(<r)-times (see [8]), we obtain 

sup£(Ao,a,u) < C^A^^I^I-^l^l-d+^l^l-d+^l^lllogr?! 2 . 
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Using that \I p \ = \I V \ > 1, \I US \ < v + 1 from [8] and that \I p \ + + |/ us | + \I CS \ = k — t, 
we obtain 

su V E(A ,a,u) < C7(A^- 4K ) dog(CT) |log7 ? | 2 

u 

if k < 1/64 and A is sufficiently small. This proves Theorem 17.11 □ 



8 Error terms: Proof of Theorem 4.1 



The main contribution to the wave function ipt in (|5.2p comes from the fully expanded 
non-recollision terms with r = 0, i.e ipfk ™" - Here we show that the contribution of all 
other terms is negligible. Our result can be summarized in the following Theorem which 
will be proven in Sections [8] and El 

Theorem 8.1 Assume (|7.ip with k < 1/9800 and a sufficiently small 5. Then, as 

mSff = o(X 4+2K+2S ) (8.1) 

for the following choices of the parameters: = rec,r = 0, 1}; {# = nest,tri,r = 1} 
or = last, r = 2}. Furthermore, for k = K and r = 0,1, 

n4?;K r f=o(\ 4+2K+2S ) (8.2) 

and for k < K 

E||V^' n l 2 =o(A 2K+25 ). (8.3) 

From this Theorem and from (15. 2p , Theorem 14.11 easily follows by using the unitarity 
estimate on the truncation of the Duhamel formula 

f* da e^-'^ips 2 <t 2 sup Us\\ 2 ■ (8.4) 

Jo 0<s<t 

Note that this estimate effectively loses a factor of t by neglecting the oscillation on the 
left hand side. (See Section 4 of [9] for more details.) 

Proof of Theorem \8.1\ We start with a resummation and symmetrization for the core 
indices that identify the non-repetitive potential labels in a sequence of collisions 7. We 
repeat the Definition 4.2 [9] here: 

Definition 8.2 (Core of a sequence) Let 7 G T n , then the set of core indices 0/7 

is defined as 

C e (7):={jG%) : Ij^lijori^j) 
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and we set 0(7) = \I^' re (^)\. The corresponding 7,- labels are called core labels. The 
subsequence of core labels form an element in T™ r , i.e. a sequence of different li d -labels. 
The elements of 

C(7):=^\[C re (7)UZ^( 7 )] 
are called non-core potential indices. 

In other words, the core indices are those skeleton indices that do not participate 
in any recollision, gate, triple collision or nest. Given our stopping rule, the number 
of non-core potential indices and ^-indices together is at most 4. The number of core 
indices c is related to the number of skeleton indices k as follows 

{k if jf = nr, last 
k - 1 if # = triple (8.5) 
k — 2 if # = nest,rec 

8.1 Resummation and symmetrization 

Let r = t(t) = (t%, ... ,r c ) G T™ r denote the core labels of the sequence 7. We rewrite 
each error term by first summing over the core labels. When computing E||^>^j^j|| 2 = 

E J tp(. . -)ip(. ■ .) by using the expansions of both ifi and ip, the core labels of ip and ip are 
exactly paired, the pairing is given by a permutation a E S c . The location of non-core 
indices within a sequence is encoded by a location code, w. The set of location codes, 
W, depends on c and r. Having specified the location of the r gates/#-indices, we 
introduce another code, h £ {g,0} r , called gate-code, to specify whether there is a gate 
or a 8 at the given location. By using a Schwarz inequality, we symmetrize for the 
location codes in the estimate of E||^|^*|| 2 , so both tp and ip have the same location 

code, w. The two gate codes, h, h' , corresponding to tp and if) are not symmetrized, since 
we still have to exploit the cancellation between the gates and #'s and this effect would 
disappear after a Schwarz inequality. See Sections 4.2 and 4.3 [9j for the details. 

For a given r, a given number of core indices c, a given permutation of core indices, 
a £ &c, a given location code w and for given gate codes, h, h! we define a partition Do 
of the joint index set I n U /„/ of ip n ip n ,. The partition Do lumps exactly those indices 
that are required to carry the same potential label by the prescribed structure. 

Some of the non-core indices may have further coincidences (e.g. two gates in the 
expansion may incidentally share the same potential label, creating a lump of four el- 
ements). This defines a new partition, D, that is the coarsening of Do, in notation: 
D y Do- Since the total number of non-core indices is bounded, the number of possible 
D is also bounded. The t?-indices remain always single. Most lumps of D have two ele- 
ments, and some gates may form quartets or sextets, their number is denoted by 04(D) 
and qq(D). Higher lumps do not appear. 
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Let D* C D denote the collection of non-single elements of D. A distinct potential 
label is selected for each element of D*, in particular the connected graph formula is 
applied for the index set D*. Let A £ »4.(D*) be a partition of the set D*. We define 
P(A, D) S V n ,n' to be the partition of I n U /„' whose lumps are given by the equivalence 
relation that two elements of I n U I n / are P(A, D)-equivalent if their D-lump(s) are 
A-equivalent. The single lumps of D remain single in P (these are the 9 indices). 

With all these notations, the following bound was proven in Proposition 4.3 of [9]: 

Proposition 8.3 Under the choices of parameters r in Theorem \8.1\ let c be given 
by (8.5\) . and let W = Wc be the set of location codes. Then 

Ell^gjll 2 < n E E E E^ (D) E c(A)V5)CP(A,D)), (8.6) 
weWaee c h,h'e{g,e} r d^d AeA(D*) 

where := E?)J are the moments of the single site random potential in (|1.4p and 

fDi 04(D) £>e(D) 
rn yy ' := mi m G . U 

The hidden parameters in the definition of the X^-value of P(A,D) (see f)6.10p ~ ()6.1ip ) 
will be chosen £ = and Q = 1 thoroughout the entire section. 



8.2 Splitting into high and low complexity regimes 



For a coarsening D >~ Do, the partition D* contains all core elements of Do- Any 
partition A G A(D*) can be naturally restricted onto the core elements and can thus be 
identified with a partition of I c (see Section 4.4 [9] for a precise definition). We denote 
this restricted partition by A. 

The sum (|8.6p will be split into two parts and estimated differently. We set a threshold 
q > 3 for the joint degree q(A, a) of a and A and obtain 



mU'Sw 2 <(/) + (//)+ o(a x / 2 ) 



with 



(/) := \W\{c+Af E E E E < (D) E bu P {|V W (P(A,D))c(A) 



w£W h,h' cree c D^D 



A'eAc 
q(A',a)>q 



(8.7) 
A = A'} 



where the supremum is over all possible A £ _4(D*) whose restriction A is the given 
partition A'; and 



iH) := \W\ E E 



E E ^ (D) E V w (P(A,D))c(A) 



Ae^l(D*) 

q(A,a)<q 
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The error term 0(X 1 ^ 2 ) comes from replacing V,°% (• • • ) with (• • • ) by using Lemma [6.51 
The high complexity regime, term (I), can be estimated by using Proposition 17.61 
By applying Operation I (Lemma 16, ip . we can break up all non-core lumps into single 
lumps, we then can remove them by Operation II (Lemma l6.2p .. The remaining partition 
is identified with P(A, a) on the core indices, I^ ore u Iff re . 

We lose at most a factor (Ar/ -1 ) 8 since Operation II is applied at most 8 times. 
Unlike in the continuous case, [9], Operation I does not cost extra factors. Following 
the argument of Section 4.5 of [9] and using the bound from Proposition 17.61 instead of 
its continuous version (Proposition 9.2. of [8]) used in [9j, we see that (I) in (|8,8p is 
negligible in the sense of Theorem 18.11 if 

•7-192 , 

k < — (8.10) 

UAq + 288 V ; 

and 5 < 5(k) is chosen sufficiently small. 

In the low complexity regime we will use the special structure given by the recollisions, 
nests, triple collisions and gates. The combinatorics in (|8.9|) was estimated in Lemma 
4.5 [9] which we recall here: 

Lemma 8.4 For any q G N, c < K and structure type we have 

sup su p Yl ' c ( A )i - (cw 3<?+3 

9(A,<t)<9 

where we recall that Dq depends on (#,c,a,w,h,h'). □ 



The size of the individual terms in (|8.9p are estimated in the following Proposition 
whose proof will be given in Section [9j This is the lattice version of Proposition 4.6 [9] 
and most of these estimates substantially differ from their continuous counterparts. The 
proof is given in Section [9j 

Proposition 8.5 We assume (I7.ip and we assume that the initial condition tpQ satisfies 
i3.20\) for some A > 0. Let a G 6 C , w G Wc , h, h! G {g, 9} r , where # and r vary in 
the different cases and let D >- Do (#, c, a, w, h, h'), 

1) Let A G *4(D*) such that q(A,a) < q. Then the following estimates hold: 

(la) [Many collisions] Let # = nr, r = 0, 1 and c = K, then 

|K(P(A,D))| <C q X 6 2 K . (8.11) 
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(lb) [Recollision] . Let # = rec, r = 0, 1, then 

|K(P(A,D))| < C 9 A^- 10K |logA|°« (8.12) 

(1c) [Triple collision] Let # = triple, r = 1, then 

|K(B(A,D))| <Cn 6 |logA|°« (8.13) 

2) Now let A' 6 A c be given. Then the following hold: 
(2a) [Non-repetition with a gate] Let = nr, r = 1, then 



sup I 

(T,W ' 



E E V(P(A,B))c(A) <C\^~ l0K - o{S) . (8.14) 



h,h'£{g,0} r D>-Do as.4(d*) 

A— A' 



(2b) [Last] Let # = last, r = 2, then 



sup 

<T,W 



E E E K(B(A,D))c(A) < CA 5 -i K |logA| W (8.15) 



h,h'£{g,0} r D^Dq A 6 .4(D*) 
A— A' 



(2c) [Nest] Let # = nest, r = 1, then 



sup 

a.w 



E E E K(B(A,D))c(A)| <CA^- 10K |logA|°« . (8.16) 

i'e{g,6} D^D AeA(D*) 

A— A' 

All constants C depend on A. 



Combining Lemma [8.41 with these estimates, and using \ W\ < K 2 , we obtain that the 
contributions of the error terms from (II) to E[|?/>^j^|| 2 (see (|8.7p ) satisfy the bounds 
in Theorem 18.11 if 

k < (8.17) 

48q + 368 V ; 

and 5 is sufficiently small. Combining this with (|8.10p and optimizing, we see that the 
biggest k that still guarantees a solution for q is around k ~ 1/9800 with q = 195. □. 



9 Proof of Proposition 18.5 
9.1 Many collisions 

The proof of Case (la) is very similar to the proof of Case (la) of Proposition 4.6 of 
[SJ; actually the estimates are sharper because Operation I does not carry an additional 
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factor A = 0(X 2dft °( <5 )) (compare Lemma 9.5 of [8J with (|6.6p ). The proof of the key 
estimate, 

I dp\K(t,p,k)\ 2 KiCaX- 2 ^- 1 

(0 < a < 1, k > TA _K_<5 , t = T\~ 2 ~ K , see Lemma 5.1 of [9]), is also simpler in the 
discrete case since no additional large-momentum cutoff factors need to be introduced. 
The only estimate that is somewhat weaker in the discrete case than its continuum 
counterpart is (|6,14p (compare with Lemma 7.1 of [8]), but for n + n' > 2K it still gives 
a negligible error. We do not repeat the details here. 



9.2 Recollision and triple collision 

The estimates for (lb) and (lc) substantially differ from their continuum counterparts 
(Proposition 4.6 of [9]). The main reason is that the two-denominator bound (|7.19p 
is weaker by a factor t/ -3 / 4- °M than its continuum version, (10.28) of [8j. Taking 
into account the A 2 -factor associated with these two denominators, each application 
of the two-denominator bound will improve the estimate by A 2 r/ _3//4_ °^ K ^ = r/ 1 / 4-0 ^ 
compared with the ladder term. The corresponding improvement in the continuum case 
is rj 1 ' ^. 

As we demonstrated in [9], a recollision Feynman graph can be evaluated by apply- 
ing the two-denominator bound twice (recall that due to symmetrization, the recollision 
graphs have two recollision; one on the ip- and one on the tp-side of the graph). The 
corresponding gain, (ij 1 ^ ^) 2 , has easily beaten the additional factor t from the trunca- 
tion (|8.4jl . In the discrete case, the recollision gain, (^V 4 " ^)) 2 , would not be sufficient. 
With the estimate (|7.19p at hand, one would need to expand at least up to five recol- 
lisions, which would mean many more terms to estimate individually. Since this route 
is quite lengthy, although certainly possible, we follow a different argument based upon 
the following 

Lemma 9.1 (i) For arbitrary a, it holds that 

f dpdqdr 
J 3 J T 3 \a — oJ(p) — ir}\ \a — Zo(q) — irj\\a — oJ(r) — vq\ \a — Zo(p — q + r — u) — vq\ 

< CA~ 3 / 4 - 3k | log??! 4 . (9.1) 

(ii) For any A > rj, there exists C\ such that for any a £ [0, 6] with |||q||| > A (recall 
the definition from \3.19\) ), we have 



sup 

u 



dpdqdr 



\a — oo(p) — irj\\a — tu(q) — in\\a — u>(r) — in\\a — uj(p — q + r — u) — in\ 

<C A \- 4K \lo gri \ u . (9.2) 
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Figure 5: Estimate of a two-sided recollision graph 




Figure 6: Estimate of a one-sided recollision graph 



Proof. The first statement (|9.ip in Lemma 19.11 is a direct consequence of (|A,24D , 
CQ5]) , (1X291) and lETR)!) . For (jn^D, we first use (b\~24j) then we use the Four Denomi- 
nator Lemma 13.41 □ 

Now we start the proof of (lb)-(lc) of Proposition 18.51 Following Section 5.2 of 
[9], Operations I, II and IV can be used to break up the partition into a trivial one 
and remove all gates and 0's. The total cost is at most A _2k | log A| 2 . Then (|8.12p and 
(|8.13p follow from the Propositions 19.21 and 19.31 exactly as the proof of Cases (lb), (lc) 
in Proposition 4.6 in [9 J followed from Propositions 5.2 and 5.3 of [9]. 

Proposition 9.2 Consider the Feynman graph on the vertex set Vk, k > 3, choose 
numbers a,b,a',b' € Ik such that b — a > 2, b' — a' > 2. Let a be a bijection between 
Ik \ {a, b} and Ik \ {a' , b'}. Let P* be the partition on the set Ik U Ik consisting of the 
lumps {j,a(j)}, j G Jfc\{o, 6} and{a,b}, {a',b'} (Fig.{5\). We assume \3.2Q\) . Then for 
the truncated version 

sup£*(P*,u) < CAx- 8K |l 0g A| 0(1) . (9.3) 

u 

We also need a "one-sided" version of this estimate (Fig. ED . 
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Proposition 9.3 Consider the Feynman graph on the vertex set Vkk-2i k > 3, choose 
numbers a,b £ Ik such that b — a > 2. Let a be a bijection between Ik \ {a, b} and Ik-2- 
Let P* be the partition on the set I^Ik-2 consisting of the lumps {j, o~(j)}, j £ Ik\{a, b} 
and{a,b}. We assume \S.2U]) . Then 

sup£(P*,u) < CA 1 -I K |logA|°W , (9.4) 

u 

and for the truncated version 

sup£*(P*,u) < CA 3 -i K |logA|°« . (9.5) 

u 

Proof of Propostion \9. 6 A We use p = (pi, . . . ,Pk+i) and their tilde-counterparts to 
denote the momenta to express 

■id 



/ia 
dad/5 S(a,/3) 
-id 



with 



l(a,/3) := X 2k j dpdpJJ 



\ l« - ^{pj) - vq\ |/3 - u(pj) + ir]\ 

X<5(Pfc+l ~ Pk+i)0~(pa+i ~Pa + (Pb+l -Pb) ~ Ua^j0~( K ~ Pa'+i + Pa' ~ (Pb'+i ~ Pb>) ~ ^a') 
k 

x H S (Pi+ 1 - Pi - (P*U)+i ~ Po-0')) ~ u i) l^o(pi)| 2 



3=1 



where the u-momenta are labelled as u = (ui, . . . , Ub-i, Ub+l, «fc 5 &> a ') 
We first prove the case k = 3, when 



/4ci r r 

dad/3 / dp 2 / 
-4d ■/ J 



dpi dp2 dp 3 

|a - cJ(pi) - irj\ |a - uJ(p 2 ) - irj\ |a - uJ(p 3 ) - irj\ 

1 



|/3 - w(pi) + irj\ |/3 - w(p 2 ) + iv\ \P ~ w fe +P3 -P2 + *4) + iv\ 
The integral over p 2 is performed by using (|7. 191) 



dp C?r 3/4_K | lo. 



< ^ ,'^ 7?l (9.6) 



2 



|/3 - w(p — u) + |/? - w(p) + irj\ 



collecting |||p3 — p 2 + Uil 1 . Then d/3, dp 2 , dp3,dpi can be performed (in this order) 
using (|A,29P to integrate the point singularity. The result is 

£*(P*,u) < CA§- 4k | log A|°^) . 
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From now on we can assume that k > 4. 

The integration will be split into four domains, depending on whether |||a||| and |||/3||| 
are smaller or bigger than A, where we recall that on the support of V'o(p) we have 
|||w(p)||| > 2A. By symmetry we can assume that |||a||| > |||/3||| and we have effectively two 
cases: 

£*(P*,u) < (7) + 2(11), 

where 

(I) := / dad/3 l(|||a||| > A, |||/?||| > A)E(a,/3) 

J~Ad 

/•Ad 

(II) := / dad/? 1(|||/3||| < A)~(a,/3) . 



Term (I): \\a\\ > A, \\/3\\ > A. 

Similarly to the proof of Proposition 5.2 in [9] , we partition the set of all p, p momenta 
into two subsets of size k + 1 each: 

A := {pi,p 2 , • • • ,Pb-i,Pb+i, ■ ..pk+x,Pv] 

B := {p!,p 2 , ■ ■ ■ ,pb>-i,Pb'+i, ■ ..pk+i,Pb} ■ 

Note that sets A and B are obtained by exchanging p^ and by i n the sets {p} and {p}. 

It is straightforward to check that all A-momenta can be uniquely expressed in terms 
of linear combinations of the 7>-momenta (plus the u-momenta) and conversely. In 
particular 

Pb-l =Pb~ (Pcr(6-1)+1 - Pa(b-1)) ~ u b 

Pb'-i = Pb' ~ (Pm+1 - Pm) + u m with 771 := a' 1 (b' - 1) . 

The letters on Fig. [5] indicate the indices of the correspoding p or p momenta. 
By the Schwarz inequality, we have 

k j 11 

n \ a - ufa) - vq\ |/3 - oj{pj) +irj\ ~ 2 ^ + ^ ^ ^ 

fa) *— TT - - TT 

J- -L \rv — 777 f m.A — in\ I R — i ,)( ri ;i\ A- in\ J- J- 



5=1 6-1,6 1° ~ U ( Pj ^ ~ ir? ' ^ ~ U ^ + ^ i=2 1° ~ U ^ ~ ^ 



^ ^ I-I \rv — (,)(rt-\ — rn\ I R — !,)( rt -i\ -A- in\ I-I 



__^_ lb \a-^(Pj)-iv\\P-^(Pj') + iv\ j± \P ~ v(pj) + iri\ 2 

jjtb'-l,b> 
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(with a little abuse of notations we used f for 1, b' — 1 and 6' when j = 1,6 — 1 and 6, 
respectively). 

Integrating out all i?-momenta in (a) and all ^4- momenta in (b), we have (with 
m:=a- x {y -I)) 

Inl<4d lfll<4ii i'=l 



|a|<4d l/3|<4d J = 



a - - irj\ |/3 - w(pi) + 

1 1 

\a - uJ(pb-i) - ir]\ IP ~ v(Pb' ~ {Pm+l ~ Pm) + u m) + 

k 



(9-8) 



n 



\a - Uj(p b+ i -p a + Pa+l ~ U a ) ~ Vq\ \P ~ 0j( Pb >) + IT]\ jj^ \at ~ U)(pj) - IT]\ 2 ' 

The integration of py is performed by using (19. 6p with v := p m +i — Pm — u m . We also 
integrate out d/3 and collect C\ log 77 1. 

Now we integrate out all pj's with j ^ l,a,a+l,b—l,b,b+l,m,m+l, each collects 
a factor A~ 2 (l + Cq\ 1 ~ 12k ) by using Lemma 13.31 Moreover, we remove the square from 
the remaining five squared denominators with j = a, o + 1, 6+ 1, m, m+ 1 at the expense 
of rj^ 1 each. Let S := {1, a, a + 1, b — 1, b + 1, m, m + 1}, then we obtain 

(/) < CA 16 rT 5 - 3/4 - K (l + OjA 1 - 12 ^- 8 ! log 7,1 / da / (J! (9-9) 

M<4d J jes 

n 



|p m+ i - p m - ii m ||| |a - w(pfe + i - Pa + Pa+i - u a ) - irj\ J-f |a - w(pj) - ir/| 

Strictly speaking, this argument is valid only if the indices 1, a, a+1, 6—1, 6, 6+1, m, m+1 
are all distinct and 6, m 7^ k. If there is some coincidence, then we may use Lemma 13.31 
for more Pj's, but save an factor each time, so the estimate still holds. In particular, 
this is the case when 3 < k < 6. 

Suppose for the moment that {m,m + 1} <f_ {a, a + 1,6 + 1}. In this case, we can 
integrate out all pj, j £ S \ {a, a + 1, 6 + 1}, using (|3.16p and (|A.29|) from the Appendix 
and the point singularity disappears at one of the integration. We obtain, by using 
k<K < A-t 1 - 12 *), 

(I) < CAf^llogr/l 4 / da l(|||a||| > A) / dp a dp a+1 d Pb+1 

I a I <4d 

n 



\a - Uj(p b+1 ~Pa+ Pa+l ~ U a ) ~ iff] \(X ~ Uj( Pj ) - irj\ 



3= 
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The proof of (|9.3p is then completed by using Lemma 19.11 

Finally, we have to discuss the case when {to, to + 1} C {a, a + 1, 6 + 1}. This can 
happen only if b — a = 2 and to = a + 1 (i.e. <y{a + 2) = <r(6) = b'). In this case the point 
singularity |||p m +i — p m — , u m || _1 would appear as an additional factor in the integrand 
with four denominators in Lemma 19.11 To avoid this situation, we choose exchange 
momenta different from Pb,Pb' when defining the sets A and B. 

If b' — a' > 2, then we can choose pb and p a '+i- The role of to will be played 
to* := a~ 1 (a' + 1). Since b' — a' > 2, we have to ^ to*, in particular {m*,m* + 1} <f_ 
{a, a + 1, b + 1} and the previous proof goes through. 

If b' — a' = 2, then we can use p a i or py+i as a tilde exchange momenta instead of 
pb', unless a' = 1 or b' = k, respectively. The role of to is played by a~ 1 (a' — 1) and 
<7 _1 (6' + 1), respectively. 

Finally, if b'—a' = 2, a' = 1, b' = k, then k = 3 and this case was already investigated. 

Term (II): \\f3\\ < A. 



From the support property of ip$ we know that 

l^obi)! 2 



\l3-uj{pi) +irj\ 



< C A (9.10) 



so the /^-denominator with p\ can be eliminated. 

We will again use a Schwarz inequality but we keep four a- and four /3-denominators 
on the first power, the corresponding index sets are A C Ik and B C Ik- 

II i -i ^ /- \ , ■ I < A-'(a) + A'(6) (9.11) 

q\ TT 1 TT 1 TT 1 

Now we explain how we choose the four-element sets A and B. 

We set A := {1,2,6 - 1,6} if a > 2. If a = 1, then we set „4 := {1,2, 3, A;} For 
the set B we set B := {1,2,6' - 1,6'} if a' > 2. If a' = 1 and 6' < k, then we set 
B := {1, 6' - 1, 6', fc}. Finally, if a' = 1,6' = fc, we set again B := {1, 2, 6' - 1, 6'}. The 
exchange momentum is always py from the tilde- variables. For the non-tilde variables 
we use pb if a > 2 and p2 if a = 1 as exchange momentum. The sets A and -B are defined 
as before: A contains all p-momenta and B contains all p momenta, except the two 
exchange momenta. 
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For simplicity, we will neglect all u momenta in the formulas below, it can be checked 
that they play no role in the arguments. 

First we compute (a), expressing everything in terms of A- momenta. The denomi- 
nator |/3 — u>(pi) + ir/l^ 1 disappears by (|9.10p . We express p2 = P\ + p c - i(i)+i — Po~i(i) 
(with the understanding that if p^ appears, it has to be reexpressed as Pb+i —p a +Pa+i) 
and we use py-i = py — (p m +i — p m ) (with m := <r~ 1 (6 / — 1)) as before. 

Suppose first that a' > 2. Then only two of the four /^-denominators contain py, so we 
can perform the py integration by ()9.6p . collecting ^- 3 / 4_K . Then the last /3-denominator 
is eliminated by the d/3-integral. For the a-denominators we proceed similarly as before 
in (19. 9D . We integrate out all but (at most) three squared denominators by Lemma 13.31 
and reduce the square to the first power in the remaining (at most) three denominators 
at the expense of r/ _1 each. Finally we use ipTIjl . The result is (a) < A 5_9k | log\\°^. 

If a' = 1, b' < k then we observe that pk is independent of py, so the same argument 
can be used as for a' > 2. 

Now we assume that a' = 1, b' = k, then both p2 and py—i depend on py. If 
I a I > A, then we estimate the p2 denominator by t?" 1 , integrate out d(3dpy at the 
expense of ^^ 1 / 2 ^ 2re using (|7.20p . After reducing the squares of the a-denominator we 
can use (|9.2p . The result is A 5_10fc | log Al ^ 1 - 1 . If |||a||| < A, then we can also remove the 
\a — uJ(pi) — zry| _1 denominator by (|9.10j) . We remove \a — u]{pb+i — p a + Pa+l) ~ ii]^ 1 
by supremum norm and note that this was the only a denominator that may have 
contained p\. If py-\ does not depend on pi, then we can integrate out dpi at the 
expense of | log A|, then we integrate dfidpy by (|7.20p and finish the argument as before 
to collect A 5 ~ 10k | log A) ^). If py-i depends on pi, then one can check that 

Pb'-i = Pv - pi + (...) 

and 

P2 =Pi - py + (■■■), 

where (. . .) refers to further A-momenta. Thus we can change integration variables, 
instead of dpidpy we consider d(pi — py)dpy. We first integrate out dpy (one j3- 
denominator is eliminated), then df3d(p\ — py) using (I7.20p . The net result of all cases 
is 

(a) < A 5 - 10K |logA|°W . (9.12) 

Now we turn to the estimate of (b). We express everything in terms of l?-momenta. 
We start with the case a > 2. Only two of the four a-denominators depend on pj,, so 
we can apply (|9.6p to perform the dpb integral to remove them. Then we estimate the 
denominator \(3 — uj(py + i —p a ' +Pa'+i) +iri\~ l by the supremum norm, thus removing the 
only /3-denominator that may depend on pi . Finally we use (|7.20p to integrate out dadpi 
if p2 depend on pi, if not, then the estimate is even better. We collect A 7 / 2 ~ 6k | log A|°W. 
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Now consider the case a = 1. We first assume b < k, then p^ is indepedent of P2 
(when expressed in terms of I?-momenta), so we can use (|9.6|) . Then we again estimate 
|/3 — uj{py + i — p a ' + Pa'+i) + ir/l" 1 by the supremum norm, thus removing the only j3- 
denominator that may depend on p\. If pk depends on pi, we use (|7.20p to integrate 
dadpi, in the other case the estimate is better. We again obtain A 7 / 2 - 6k | log A|°W. 

Finally, we consider the case a = 1, b = k. We estimate \/3 — oj{p~y + i — p a i + 
Pa'+i) by supremum norm so no j3 denominator can depend on pk+i- We express 

Pk = Pk+i —&1+P2- If P3 depends onpk+i, then we perform dpk+i using (j9.6|) . then we 
can perform dp2 removing one a-denominator and finally we integrate out a to remove 
the last a-denominators (with pi). The remaining /3-denominators are independent and 
we collect A 9 / 2 ~ 4k | log A) * 1 ). If p% does not depend on pk+i, then we integrate out dpk+i 
to remove one a-denominator (with p^) and collecting | log A| . Then we integrate p2 using 
(|9.6p and finally perform the a-integration. The result is again A 9//2_4k | logX\°^. 

In summary, we obtain 

(b) < A 7/2 - 6K |logA|° (1) , (9.13) 

and together with (|9.12p and optimizing for p in (19. lip we obtain (|9.3p . This completes 
the proof of Proposition 19.21 □ 

Proof of Proposition \9.SX This proof is similar to the previous one but simpler. We 
choose the set of A and B momenta are as follows: 

A := {p!,p 2 , ■ ■ ■ ,Pb-i,Pb+i, ■ ■ ■ ,Pk+i} , B := {pi,p2, ■ ■ ■ ,Pk-l,Pb} 
If I a I > A, then the Schwarz estimate is the following 
fe+i fc-i 1 

TT 1 =7^ TT T-o ~7 < A » + X ~ P ( b ) 

11 U _ —(nA — vn\ 11 \R — Mfn.l 4- tn\ w v ' 



ufa) - irj\ AJL \p - ufa) + irj\ 



1 1 fe+1 1 

/ \ A TT A TT 1 

I R — tn(ri-i \ — in\ 11 I rv — T7i(r>A — in\ 11 



|/3 - uj(pi) - i<q\ *■ |a - w(p 3 -) - m?| ^ |a - u(pj) - irj\ 2 



1 ^ 1 



^ ' I/ 3 - u(pi) - H JJ 6 | a - w(pj) - ir/| n _ ^(p,.) + ^|2 ■ 

To estimate the integral of (a), we express p\ = p\ and p& = pb+i — p a + Pa+i — ""a, so 
every term in (a) will depend only on j4-momenta. We first integrate d/3, then integrate 
all pj, j ^ a, a + 1,6 + 1 and reduce the square of the remaining denominators to the 
first power. Finally we use Lemma Ell The result is CA p ~ 7k | log A|° (1) . The integral 
of (b) is even easier, after expressing p% = pi, we can integrate it out in any order 
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with an estimate CX 2 p \log X\°^ . After optimizing for p this gives 

A l-W2|l 0gA |0(l) 

as announced in (|9.4p . 

If I a I < A, then we use 

fe+l k-i 1 

1 = 1 |a - UJ{ Pj ) - i V \ 1 = 1 |/3 - wfo) -h 

1 t-t 1 1 



: n \R-,,,(aA-i.'n\ n u _ <r,( nj \ _ 7 -«i n 



f 2 \P ~ u(Pj) ~ H , =1 b i b \ a - ^iPj) - H |a - w(pj) - ^| 2 

1 1 fc_1 1 

(6) IJ 2 1/3 - w(pj) - i= JJ 1)6 la-wCpj)-^! 2 |/3-w(^)+^| 2 ' 

We use (I9.10|) to eliminate \a — u(pi) — ir/| _1 . 

In the term (a) we express everything in terms of A- momenta. We estimate \a — 
Za(pb+i — p a + Pa+i) — iri\~ l by rj~ l . Thus no more a-denominator depends on p\. 
Depending on whether p2 depends on p\ or not, we can use (|7.20p or subsequent dpi 
and d/3 integrations to remove all /^-denominators. The remaining a-denominators are 
independent and we collect 

{a)<\- x -^\\og\\°^ 

In the term (b) we use the .B-momenta for integration. Here only two a-denominators 
depend on p , so we can use (I7.20P to perform dp da, then all /3-denominators are 
integrated independently. We obtain 

(6) < A 3 ~ 2k | log A|° (1) . 

After optimizing for p, we obtain a bound smaller than (|9.4p . 

The proof of (19. 5j) is the same, but the last squared denominators are missing, this 
is where the gain A 2 comes from. □ 



9.3 Cancellation with a gate 

The proof of (2a)— (2c) of Proposition 18.51 depends on a cancellation mechanism between 
a gate and a 8 label. More precisely, if two Feynman graphs differ only by replacing a 
gate with a #-label, then their sum is by a factor A 2 ?? -1 / 2 smaller than the E- value of 
the two partitions individually. Moreover, this cancellation effect is local in the graph: 
if another gate/9 pair occurs somewhere else in these graphs, doubling their number, 
then the sum of these four Feynman diagrams is smaller by a factor (A 2 ?^ 1 / 2 ) 2 . For 
the general statement, see Lemma 5.5 in [9]. Although this lemma is formulated for the 
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continuum model, taking the collision function B = 1 and considering all momemtum 
integrals in T 3 , the proof of Lemma 5.5 goes through for the lattice case as well. 

The detailed proofs of (2a)-(2c) follow the arguments of Section 5.3.2-5.3.4 of [9] line 
by line and will not be repeated here. We only point out the three minor differences: 

(i) The A-exponent in the estimate on Eu\(A, a, u) given in Corollary 17.51 differs from 
its continuum counterpart (9.4) of [8]; 

(ii) In the continuum model, each application of Operation I costs a factor x~ 2dK -°(°) 
(denoted by A in 012]), this loss is absent here; 

(iii) The A-power in the estimates ()9.3|) — f)9.5|) are weaker than their continuum ana- 
logues (Propositions 5.2 and 5.3 of [9]). 

These changes account for the somewhat different A-powers in (2a)— (2c) of Proposi- 
tion [83] compared with Proposition 4.6 of [9j. □ 



10 The main term: Proof of Theorem 4.3 



For simplicity, all results in this section are written for d = 3, the calculation for higher 
dimensions is similar. We follow a different path than in the proof of the analogous 
theorem in Section 6 of [9]. Due to the uniformity of the Boltzmann collision kernel 
(|2.19p . we can circumvent the reference to the Boltzmann process and we identify the 
heat equation by a direct computation. 

As in Section 6 of 9]. we start with the identity 

W x (t,k,0)= [ V£(A ,6(£, ■))<!£, k>l 

J{2T/e) d 

with Ao being the trivial partition on 1^, where we chose the function Q(v) in the 
definition of V° to be ^-dependent, namely Q(v) = Q((v) := 0(£,v) (see (|6.10p and 
(|6.13p for definitions). 

Analogously to the argument in Section 6 of [9], the d£ integration can be restricted 
to the regime {|£| < \~ s } with a negligible error (even after summation over k): 

W x (t,k,0)= + E °k--= f ^(A o ,0(£,-))d£, (10.1) 

l<k<K l<k<K J 

where we used the notation 

(*{■■■)<%:= I (•••)l(l£|<A- 5 )d£. 

J J(2T/e) d 
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The bound ([UTT]) follows from LemmaESI from |^(A Q , •)) | < \\Qt\\oo sup €jU E{a = 
id, u), from the uniform bound (|7,10p - (|7.12p on E[a = id, u) and from the fast decay of 
||Q e || 00 = sup 1 ,|0(e,u)| in£. 



Writing out the definition of H£ more explicitly, we have 

, fc+i 



?o _\2k 
J fc — A 



r- 

fc+l j— 



II + f)^(/?,^ - y)] • (10.2) 

The estimates of the error terms were performed with the choice r\ = A 2+K . However, 
H?, given by (|10.2p . is independent of r\. Therefore we can change the value of r] to 
77 := A 2+4k for the rest of this calculation and we define 

R(a, v) := R v (a, v) , with 77 := A 2+4k . 

We recall that the restriction of the dad/3 integration in (|10.2p to any set that contains 
{a, (5 : \a\,\(5\ < 4d} results in negligible errors, even after the summation over k 
(LemmaESJ). We will consider the set D := {(a,/3) : \a + f3\ < 8d, \a - (3\ < 8d}. We 
denote by the version of z? k given by formula (|10.2p with the dad/3 integrals restricted 
to D, 

Sfc .= A 

then 



f f 7— Integrand from (110. 2p 
J Jd ( 27 L 

53 \Z° k -E k \ = o(l). 



Kk<K 



We also remind the reader that this argument does not apply literally to the trivial k = 
case, when the da d/3 integral in (|10.2p gives free evolutions and this term is computed 
directly: 

H : = J* d£dv e it£ ^< v) e 2tx2lme ^ v)W^(e£, v) + o(l) . (10.3) 

The error term comes from the error term in uJ{v + — u(v — e£/2) = el; ■ Ve(w) + 

2i\ 2 l{v) + 0(A 2 e£) + 0{e 2 i 2 ). By using t\ 2 -> 00, the bound t$3J)§, and the decay of 
the observable, one easily obtains that |Eo| = o(l). 

We start with a crucial technical lemma which is proven in the Appendix. 
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Lemma 10.1 Let k < 1/18 and set 7 := (a + (3)/2. Let n satisfy A 2+4k < n < A 2+K . 
JTien /or |r| < A 2+K//4 and any / 6 C 1 (T £i ) we /iaue, 



A 2 /(P) 



(a — u>(p — r) — in) (/? — oj{p + r) + in) 
2 I' fiP) He(p) -7) 



dp 



(10.4) 



(«-/?) + 2(Ve)(p) • r - 2i[A 2 X( 7 ) + n] 



d P + 0(X 



1/2-9/t 



C 1 



We now apply this lemma to compute Introduce new variables as a := (a + (3)/2 
and b := X~ 2 (a - (3). Then 



A 2 J duT(£, v) R(a,v + y) v ~ "4] 

6 + £ A- 2 (Ve)(,).e-2,[J(a) + A^] d " + ° (A )l|T|kl ' 



(10.5) 



We now replace the product of k + 1 factors in the restricted version of (|10.2p one by 
one. We need a A 2 factor for each application of f)10.5j) . The (k + l) st factor A 2 comes 
from the change of variables dad/3 = A 2 dad6. We also define D* := {(a, b) : \a\ < 
Ad, \b\ < 8dA -2 } as the domain D in the new variables. 

Introduce the notation 



Fi(£,v) := W^fov), F k+1 (£,v) := O&v) 

and Fj(£,v) := 1 for j = 2, . . . k. Using (|10.5p with T = Fj and using that ||C||ci and 
HW^ollc 1 are bounded, we obtain by a telescopic summation that 



k<K k<K 



d£ 



D" 



dado 



fk+1 

- n 



Jt\ 2 b+2trj 



-2mFj(^,Vj) 5(e(vj) — a) 



i ., 6 + eA- 2 (Ve)(v i ) • £ - 2t[X(a) + A 4K p ^ 



(10.6) 



with 

^ M := 0(\ 1 ' 2 -^) 



k+l 



D* 



dadfe 



n 



27ri o~( e ( u j) — a) 



6 + eA- 2 (Ve)(uj) • £ - 2i[I(a) + A 4k ] 



n 



(10.7) 
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The first factor in (|10.7p is zero if a (0, 2d), otherwise it can be estimated by 
2m 8(e(v) — a) 



I 



b + e\~ 2 (Ve)(v) ■ £ - 2i[l{a) + A 4k ] 



dv < —— f vr 5(e(v) - a)dv = 1 (lOi 
1(a) J 



by using f)3. 13j) . The second factor in (|10.7p can by bounded by the Schwarz inequality 
and by (pTTB]) : 

A 2 sup / \R(a,v)\ 2 dv < 1 + 0(A 1 ~ 12k ) . 

a J 

Thus the right hand side of (|10,6p vanishes in the limit A — > since the double summation 
yields only a factor K 2 = 0(A" 2k - 2<5 ). 

Now we concentrate on the main term, i.e. on the sum on the left hand side of (|10.6p . 
We first remove A 4k from the denominator in the main term in (|10.6[) . To estimate this 
replacement error, we first notice that the large regimes of a or 6 are harmless. Due 
to S(e(vj) — a), the integrand is zero unless a £ (0, 2d). In the large b regime, each 
denominator can be estimated by (|&| + 1) _1 since \e\~ 2 (Ve)(v j) ■ £| <C 1 and there are 
at least two denominators (k > 1), so the tail of the d& integration is small. 

The integration in a is divided into two cases. In the regime where 1(a) > A 3k/2 , the 
error of the replacement in one denominator is bounded by 



C 1 1 1 1 oo 



— 5{e ^ - Q ) d, < 0(X^ 2 ) 



by using a resolvent expansion. For the complement regime, 1(a) < A 3k / 2 , we will use 
the trivial estimate (|1U.8|) . Therefore we obtain 



-2mFj(j£,Vj) 5(e(vj) — a) 

dvj 



b + eA- 2 (Ve)(uj) • f - 2i[l(a) + A K / 4 ] 

[ - 27H Fj(t,Vj)S(e(v3)-a) , . 0(x 5k/2, , nh) , (T(fA < a 3«/2n ( , n q) 
-J b + e\-*(Ve)(v J )-ti-2il(a) dv ^ +0 ^ X 0(1)1(1 (a) < A ) . (10.9) 

In the regime 1(a) > A 3k / 2 , similarly to the telescopic estimate leading to (|10.6j) . 
we can remove the A K//4 terms from each denominator on the left hand side of (|10.6[) 
one by one, the error is controlled by CK 2 \~ 5d \^ K / 2 — > if 5 is sufficiently small. For 
1(a) < A 3k / 2 we use the trivial estimate (|1Q.8|) for each denominator, and the fact that 

da l(l(a) < A 3k / 2 ) J |T^(e^t^)| <5(e(«i) - a)dv x 

< C\\Ml j 1(a) l(l(a) < A 3K / 2 )da = CA 3 «/ 2 ||^ ||L • (10-10) 
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Together with the summation and the d£-integration, this yields an error term of order 
KX 3K/2 X -8d q_ We therefore obtain 

E~ fc = o(l) (10.11) 



k<K 



, V- [* ac [ f dadb Jt\ 2 b+2t n ( TT f -ImFj&Vj) 5(e( Vj ) -a) \ 

Notice that we removed the constraint 1(a) > A 3k / 2 after having used it and we ex- 
tended the dadfe integration from D* to R 2 . The corresponding errors are negligible by 
arguments similar to the previous ones. 

Next we perform a change of variable b — > fo(2Z(a))~ 1 so that 

d6 e ux2b+2t ^ - 21(a) db e ^(a)tb+2t v 

Since trj = A 3k — > 0, we shall drop the e 2tr) factor. Denote H(v) := ■ Introduce the 
probability measure d^ a (t>) on the level surface S a = {e(v) = a} by 

h(v)dfi a (v) := (h) a = I h(v)5(e(v) - a)dv . 

Then 

f -2mF J (C,v)5(e(v)-a) f -iF&v) 

with the new variable b on the right hand side. 

From (|3,14[) . H is bounded above, i.e., sup„ < C. With this bound, we expand 

the fraction up to second order in eA _2 |£| < A K / 2 ~ 5 <C 1 for 2 < j < k 



/ 



b + e\- 2 H(v) -i-i 



d/J, a (v) 



-i 



1 — + — (b—iT 2 — + °{ {£X 



dfi a (v) (10.12) 



Since e = A 2+K / 2 , |£| < \~ s and K = 0(A~ K ~ 5 ), the last error term, 0[(e\- 2 \i\f 
even after summation in k, is negligible: 

^(eA- 2 |e|) 3 = o(l). (10.13) 

By symmetry, H(v) = —H(—v). Therefore the linear term in £ on the right hand side 
of (|10.12p vanishes after the dfi a integration since dfi a (v) = dfj, a (—v). 
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For j = 1, k + 1 we will use the following simple estimate 



I 



b + e\- 2 H{v)-i-i 

We also define the matrix 
21(a) 



b-i 



F j (tv)d^ a (v)+o(l) 



(10.14) 



D(a) :-- 



(27T) 



d/i» H(v)®H(v) 



1 



2vr$ a 



■( sin(27Tf ) (g) sin(27n;)) a 



using 1(e) = 7r3>(e). The quadratic form of D(a) is denoted by (£, Z?(a)£), £ G 

By applying (|10. 14|) to j = 1, k + 1 and ()1Q. 12[) to the rest of the du,- integrals in 
(jlQ.lljl . we obtain 

Y, W x (t,k,0)=-^ f d£ ! da21(a) / W^(e£, «i)d// («i) 
x / 0(M<K(«) / d6e 2iA2x ^ fe 



i<fe<^ 



E 



6-i 



fc+1 



1 + 



(2^)VA" 4 (£ •£(«)£) 1 



21(a) 



(b-if 



(10.15) 

fe-i 



modulo negligible errors. We shall also change the last exponent from k — 1 to k + 1 
to simplify the computation. Since e 2 A~ 4 — > 0, this modification causes only negligible 

errors. 

Since e is small, we can replace W^ (e^, v±) by W^ (0, «i) G R with a negligible error, 
by recalling that ipo S C 1 (Section 13. 2|) . in particular, 



W^ (e(,vi)dn a (vi) = g(a) + 0(e) , g(a) := [|^ | 



Let A := 2\ 2 l(a) ^ and B := 2^eA~ 2 ( 



1/2 



Suppose that we can extend the summation in A; to infinity in fjlQ. 15|) . Then 

fc+l 



E 



fc=0 

and we compute 



b-i 



k+l 



1 + 



B 1 



(b-i) 2 
(-i) I db e ltAb - 



Note that B < 1 uniformly in a. 
in CELL 

(b-i) 2 + B 2 



(b - if + i(b - i) 2 + iB 2 
(b-i) 2 + B 2 



i) 3 + i(b - i) 2 + iB 2 ' 
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The roots of the denominator are &i, &2 ; ^3 and using that fi < 1, we obtain 

bi = iB 2 + 0(B 4 ), b 2 =i + iB + 0(B 2 ), b 3 = i - iB + 0(B 2 ) . 

The roots 62, 63 give exponentially small contributions since e ltAb2 ' 3 < e~ tA = exp(— A~ re <&(e 
The main contribution comes from b\, so 



H) / dfte^' 
We have 



(b-i) 2 + B 2 



27te- tAB \l + 0(B)) + o(l) 



i) 3 + i(b - i) 2 + iB 2 



tAB 2 = (2ir) 2 e 2 \- 4 - K T{£,D(a)0 = (2vr) 2 T(£, D{a)£). 

Note that to get a nontrivial limit, the space scale e has to be chosen as e = \~ 2 ~ K / 2 
here. Thus we obtain 



Km Yl W x (t,k,0) 



(10.16) 



Kk<K 



f d£ f p21(a)( I O(e, V )d / u a ( W )) 5 (G)exp(-(2vr) 2 r(e, J D(a)0 

Since f(T,X,a) satisfies the heat equation (12.150 with initial condition f(0,X,a) := 
5(X)g(a), its Fourier transform in X is given by f(T, £, a) = g(a) exp(— (2n) 2 T(^, D(a)^)). 
Thus 

lim V W A (t,fc,0) = (2^)- X / d£ I da21(a) I d(^,v)d f i a (v)f(T^,a) 

d£ [ dvd(£,v)f(T,£,e(v)) , 



where we used the definiton of dfi a . This proves Theorem 12.11 

It remains to prove the contribution from k > K is negligible. By the residue theorem 



/ 



dbe 



itAb 



b-i 



with 



It is easy to check that 



1 + 



a 



B 2 



fc+i 



(2k) 



I, fe+1 ! 



{At) p _ tA 
k\ 



(BAtf 



k 



;io.i7) 



k,t 



kl(k + l)lk £ 



(k + 2£)l(k + 1 - £)V 
Cfc.Q = 1) Chg < 1 • 



55 



If k > K = 0(\- K - 5 ), then At/k = 0(\ s ) and (BAtfjk < X s (using that 1(a) and 
D(a) are bounded), so (|10. 17|) is smaller than 

(At) k _ tA ^kxkS 

~TT e - A ' 

which is negligible even after summing up for all k > K. □ 



A Estimates on Propagators 



A.l Proof of Lemma |3J3 



Let Hffllc-i = 1 1 (7 1 1 oo + ||V<7||oo denote the C 1 norm of a function on the torus T d . We 
start with the following Lemma. 



Lemma A.l For any a, a' real and e > e' > 0, we have in d > 3, 
1 1 



a — e(p) + is a' — e(p) + is' 
In particular, with g = 1 we have 

\e £ (a) - Q £ ,(a')\ < C 



g{p)dp < C 



e -l/2| a _ c /| +e l/2 



Iffllci ■ (A.l) 



(A.2) 



The functions $ (from (|3.1(jp ) and Q are Holder continuous with exponent i for d > 3: 
Me) - ${e')\, |0(e) -9(e')| < C\e-e'\ 1/2 . (A.3) 

We also have 

' ' ^dp^CA^IlogAl \\g\\ c i (A.4) 



a — uj(p) + irj a — e(p) + ir\ 
for any n satisfying A 2+4k < n < A 2+re . 

Proof. We shall consider the case d = 3 only, for d > 3 the proof is similar. We first 
describe a local coordinate system. Let I J + = {pj : \pj\ < 2ir/3} and P_ = {pj : \pj\ > 
7r/3}. Let g G {(±,±,±)}. We can find a smooth partition of unity \a i n the sense that 



(A.5) 



and Xa is supported in l a = n|=i ^ • Each I CT contains exactly one critical point of 
e(p). 
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Clearly, g(p) = ^2 a Xu(p)g(p) and we can prove the lemma for each a fixed with 
g replaced by gXo- We now consider the case o = (+,+,—) which corresponds to a 
hyperbolic critical point. Define 

Uj = V2sm(irpj), j = 1,2; u 3 = V2 cos(irp 3 ) (A. 6) 

We have 

dpidp2dp3 = J(u)duidu2dii3 (A. 7) 

where J(u) is a smooth, non- vanishing function on I a . In terms of u, the dispersion 
relation in I a can be written in the canonical form 

e(p) = u\ + u\ - u\ + 2 (A.8) 

Thus u = is a hyperbolic critical point. For all a fixed, we can perform a similar change 
of variables. There are two elliptic points and six hyperbolic points. The following 
argument for a = (+, +, — ) applies to all hyperbolic points. For the elliptic points, after 
the change of variables into the canonical form, the same result was proved in Lemma 
3.10 of 0. 

Returning to the case a = (+, +, — ), we need to estimate 



/ 



1 



a — (uf + «2 — ii§) + ie a' — (u{ + — uQ + is' 



f(u)du (A.9) 



where we have replaced a — 2 by a and a' — 2 by a' . Here f(u) = g(p(u))X(r(p(u))J(u), 
where the regular function p(u) is given by the inverse of the change of variables formula 

(ED. 

Define U\ = Q 1 / 2 cos <p,U2 = Q 1//2 sin</>, 

f-2n 

f*(Q,u 3 ):= /(Q 1/2 cos0,Q 1 / 2 sin^,u 3 )d0 (A.10) 
Jo 

Y(z,u 3 ) := f r{Q ^ ] dQ (A.ll) 

Let z = a + n 2 + ie and z' = a 1 + n 2 + ie'. For any fixed U3, we can follow the proof of 
(3.69) in the appendix of [6] to have 

\Y(z, u 3 ) - Y(z>, u 3 )\ < C\\f*\U log z - log z'\ + C\z - A^WSW^ ■ (A.12) 

We point out that there is a typo on the right side of equation (A. 7) of [Bj. The correct 
bound should be const. |Im ^l -1 ^ 2 !!/!!^ , which leads to (|A.12|) . 
Since 

\z — z'\ < \a — a'\ + e , 
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the last term on the right hand side of (|A. 12() . after integration in U3, is of the form 
stated in (|A.ip , For the first term, we use 

/ du 3 | log(a+ul+ie)-log(a'+ul+ie')\ = / — -^-1 \og(a+Q+ie)-\og{a +Q+ie')\ ■ 
Jo Jo Q ' 

The last integral is very similar to the equation following (A. 6) in |6]. We can follow the 
proof in [6] after (A. 6) to estimate this integral by C\z — z'\e~ 1 / 2 . We thus bound both 
terms in (|A.12|) and prove (|A.1|) . Equation (|A.3|) follows from optimizing e in (jA.ip . 
To prove (IA.4j) . we rewrite it as 



g(p)dp 



+ 



a — lo(p) + ir] a — e(p) + ir] 

A 2 (8(e(p))-e(a)) 

(a — e(p) — \ 2 @(e(p)) + irj)(a — e(p) — A 2 G(a) + ir]) 



g{p)dp 



(A.13) 



/ 



a — e(p) — A 2 0(a) + ir] a — e(p) + ir] 



recalling that uj(p) = e(p) + A 2 0(e(p)). From the Holder continuity estimate (|A.3|) . we 
can bound the first term on the right hand side by 



C\\g\\ 



\e(p) — a] 1 / 2 



A 2 



(a — uj(p) + ir]) (a — e(p) — A 2 0(a) + irj) 



dp . 



(A.14) 



Since 



\e{p) - a\ 1/2 < \u{p) - a\ 1/2 + O(A) , 
the integral in (|A,14j) is bounded by 

dp 



\a — e{p) — A 2 0(a) + ir]\ 
By the co-area formula (|2.5p . the last integral is bounded by 



dp 



\a — e(p) — A 2 ©(a) + ir]\ 



2d 



ds 



a 



X 2 Q(a)+ir]\ 



< C\logr]\ . (A.15) 



Here we have used that is bounded according to (|3.1ip . This bounds the first term in 
(|A.13p . The second term is bounded O(A) from (|A,3p . This proves Lemma lA.ll □ 



Proof of Lemma The proof for (|3.16|) and (|3.17|) is similar to the argument 

from ()A.13P to (|A.15P and we shall not repeat it here. We now prove the more accurate 
estimate ()3. 18j) : it is is similar to the proof of (2.9) in [9] (Appendix B.l). 
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From the Schwarz inequality and a change of variables, it suffices to prove only the 
first estimate of ()3.18p . We can assume that \a\ < 4d, since otherwise there is no 
singularity at all and the Lemma trivially holds. 

Recall w(p) := e(p) + X 2 9(p), 9(u) = 9(e(u)) and 9(e) = 11(e) - il(e) with 1(e) > 0. 
We have 

A 2 A 2 



\a — u(u) — irj\ 2 X 2 l(e(u)) + rj 



Im 



From the resolvent identity and with the notations e 
equal to (I) + (II) + (III), where 

rn _ A 2 1 
1 ' " A 2 J(5) + rj m a - e - i(X 2 l(a) + rj) 
A 2 X 2 (1(e) -1(a)) 



a - e(u) - \ 2 TZ(e(u)) - i(\ 2 l(e(u)) + rj) ' 

e(u), a = a — \ 2 lZ(a), this is 



(A.16) 



(II) :- 
(III) 



\ 2 l(a) + t] \ 2 l(e) + r] 



Im- 



Q: 



A 2 



A 2 T(5) + 7? 



Im 



A 2 ©(e) — ir] 
A 2 (B(a) - 9(e)) 



a 



i(A 2 T(a) + ??) a 



A 2 9(e) — ir] 



We will estimate J du[(I) + (-f-f) + (^^^)]- The main term will be the first one. In this 
term we first use the definition (|3.2|) and the estimate (|A.2|) with e := X 2 l(a)+rj = 0(X 2 ) 
to obtain 

J(2) + 0(A)1 < 1 + 0(A 3 ??- 1 ) . 



/ 



(/) du < 



X 2 l(a) + r/ 
In the second term of (|A.16|) we use (1A.3P 



\(II)\ <CX 



2/A 2 X 2 



a 



1 1/2 



e + A 2 (7e(5) -n(e))\ 2 + r] 2 ' 



where we also used that A 2 (7£(a) - 11(a)) = 0(A 6 ) < r/ 2 based upon (|A.3|) . 

To perform the du integration (recall e = e(u)), we distinguish two regimes depending 
on whether \a — e(u)\ is bigger or smaller than KqX 4 ^ for a sufficiently large fixed Kq. 
When \a - e(u)\ > K X 4 , then X 2 \TZ(a) - 1l(e(u))\ < \\a - e(u)\. Hence 



|(//)|<CA 2 



A 2 \2 



,-1/2 



\a — e(u)\ + rj 



and the corresponding integral is of order 0(A 6 77 _5 / 2 | log r?| ) using the co-area formula 
and the boundedness of <£. When \a — e(u)\ < KqX 4 , then we can trivially estimate 
|(7Y)| < C(A 2 r/ -1 ) 4 and the volume factor is given by 

/ l(|5-s| < K X A )$(s)ds = 0(X 4 ) . 
Jo 
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Therefore the integral is of order 0((A 3 ?7 x ) 4 ). 
Finally, the last term in (|A.16|) is estimated as 



|(7/I)| <CA 



. 2 /A^ 1 la-el 1 ^ 



rj ) \a — e| + 77 \a — e + X 2 1Z(e)\ + r/ 
In the regime where \a — e(u)\ > -KoA 2 (with some large Kq) we obtain 

\(iii)\ < cxH-) „ la ~f\ 9 < caV 1/2 ^ 1 



77 / (\a — e\+rf) 2 ' \ r] J \a — e\ + i] 

and after integration we collect 0(A 4 r7 _3 / 2 | log??!). ^ n the regime where \a— e(u)\ < Kq\ 2 
we have |(/iT)| < 0(A 5 ?7 _3 ) and the volume factor is 0(A 2 ). Therefore the integral is 
0(A ry -3 ). Collecting the error terms, we arrive at the proof of Lemma 13.31 □ 



A. 2 Proof of Lemma 110.11 

We can assume that / is a real function. We can also assume that \a — (3\ < A, otherwise 
the two singularities are separated by O(A) and at least one of the denominator can be 
estimated by 0(X~ 1 ) and the other one integrated out by (|3.16l) to give 0(X\ logA|). 

We use the partition of unity 1 = ^ CT Xo- introduced in (|A.5P and let f a := fxa- 
Fix a a around a hyperbolic critical point; the proof for the elliptic point is similar but 
easier. We choose the case a = (+,+,—) as in the proof of Lemma IA.1I Recall the 
new coordinate system on I a from (|A.6|) . where \uj\ < \f2 — cq with some positive cq. 
In terms of u, the dispersion relation in I a can be written in the canonical form (|A.8D . 
Introduce the notations 

Gj(u) := 2tt^2-u] Uj , j = 1, 2; G 3 (u) := -2vr^2-n 2 u 3 , 

then the vector-field G is regular on the integration domain. Note that 

(Ve)(p(«)) = G(«). 

Let 

f(u) = f a (v{u))J{u) , (A.17) 

where J(u) is the Jacobian (IA.7p . We expand in r as e(u ± r) = e{u) ± Ve(n) • r + 0(r 2 ) 
and we also use the Holder continuity of given in (|A,3D : 

e(e(p(u) ± r)) = 6( 7 ) + 0(| 7 - e(p{u))\ 1 / 2 ) + 0(\r\ 1 ' 2 ) . 
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Thus we can write the integral (|10.4p around the selected critical point as 



duf(u) 



(A.18) 



A 



a - (u\ + u 2 - ul) - G{u) ■ r + 0(r 2 + A 2 |rf/2) _ A26( 7 ) + 0(A 2 | 7 - e|V2) _ ir] 

A 

p- (ul + uj -ul) + G(u) -r + 0(r 2 + A 2 |rf/2) - A 2 9( 7 ) + 0(A 2 | 7 - e|V2) + ir] 

where e = e(p(u)), a := a — 2, /?:=/? — 2. 

The error terms can be removed from the denominators by an argument similar to 
the proof of Lemma 13.31 in Appendix IA.1I Using the resolvent identity for the first 
denominator in (|A.18P 

1 1 B 



A + B-irj A-irj (A + B — irj)(A — irj) ' 

with A := a - (ul + u\ - u 2 ) - G(u) ■ r - A 2 8( 7 ), B := 0(r 2 + A 2 |rf/ 2 ) + 0(A 2 | 7 - e] 1 / 2 ), 
and | 7 — e] 1 / 2 < \a — e\ 1 / 2 + 0(A 1 / 2 ) < l^j 1 / 2 + 0(\ l l 2 ), we can estimate the error term 
by 



B 



(A + B — ir])(A - irj) 



< 0(t?- 2 [A 5 / 2 + r 2 + A 2 ^! 1 / 2 ] ) + 0(AV 3 / 2 ) = 0(A- 3 / 2 - 8K ) . 



Then we can integrate out the denominator with (3 to collect an extra | log A| factor. 
This term, together with the A 2 in the numerators of (IA.18j) . is included in the error 
term in Lemma 110. 11 

Now we compute the main term 



M :-- 



duf(u) 



A 



a — (u 2 + u\ — u 2 ) — G(u) ■ r — A 2 G( 7 ) — ii] 

A 



(A.19) 



(3 -(u\ + u 2 2 - uf) + G(u) ■ r - A 2 6( 7 ) + ir] 
We introduce spherical coordinates in the (u\, U2) plane. Let Q := u\ + u\ and 
u\ = Q 1 ! 2 cos <j), U2 = Q 1 ^ 2 sin <j) . 
We shall use Q,4>,us = y as coordinates instead of m, U2, U3 then 

duidn 2 du3 = -dQdydcj) . 
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Thus 



M 



dydQd<t>f(Q,y,, 



X 



a + y 2 -Q- G(Q, <p,y)-r- A 2 6( 7 ) - irj 
A 

< p + y 2 -Q + G(Q, ^y)-r- A 2 9( 7 ) + iq 



(A.20) 



The integration domains are Q G [0,4], y G [— \/2, a/2], <A G [0, 27r], but since / is 
compactly supported on this domain, we can assume that Q G [0, oo) and y G E. Set 
p := 7 + y 2 , 7 = 7 — 2, and distinguish three cases (subdomains of integration). 

Case 1: \p\ < 3A. After a Schwarz inequality, estimating / by its maximum and a 
change of variable Q — > Q := Q ± G(Q, 4>,y) ■ r we estimate 



dydQd0 



A 2 -1(|7 + ^ 2 |<3A) 
5 + y 2 — C/ — A 2 G(7) — Z77| ; 



< CAV 1 f dy l(|7+y 2 | < 3A) < 0(A 



1/2-4k^ 



Case p < — 3A. The contribution of this regime, after a Schwarz inequality and 
the estimate 1(7 + y 2 < — 3A) < 1(5 + y 2 < —2 A), is bounded by 



A / dyd^ 



A 2 



dQ 



1(5 + y 2 < -2A) 



<CA 2 



dyd^dC/ 



Q 



a + y 2 — irj\ 

(A.21) 



|5 + y 2 -g-A 2 6(7)-i?7| 2 

<0(AV 1 / 2 |i og A| 2 )<0(A 1 / 2 ). 

The integration domain of Q comes from the definition of Q, Q > and r -C A 2 . 

Case 5: p > 3A. We replace G(Q,(j),y) by G(7 + y 2 ,0,y) in fA20l) . Since G is 
regular on the integration domain, the error of this replacement in M is estimated by 
0(A 1 / 2-8k ) by a similar argument as we removed the error terms in (|A.18|) . Therefore 
the contribution of the regime p > 3A to (|A.20p is (up to negligible errors) 



M 3 := 



2 h 



7+y 2 >3A 



dyd<$> \ dQf(Q,y, 



X 



a + y 2 - Q + G(j + y 2 , 4>,y) -r — A 2 0(7) - ir\ 

X 

P + V 2 - Q + G(j + y 2 , 0, y) ■ r - A 2 6( 7 ) + irj 
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1(7 + V 2 > 3A)dyd</> 



1 



a -13 + 2G(7 + y 2 ,cf>,y) • r - 2iA 2 X(7) - 2ir) 

1 



(A.22) 



AQf(Q,y, 



1 



/3 - Q + - A 2 9( 7 ) + ir] a - Q + R - A 2 9( 7 ) - irj 
with i? := y 2 + G(^y + y 2 , cp,y) ■ r and recalling that G = TZ — il. 
Lemma A. 2 Let F be a C l -function on M with compact support and let 

HQ) 



Y(z) 



o z - Q 



dQ 



for any z = a + iA with A > 0. Then 

\Y(z)-Y(z')\ < C\z- z'\\\ogA\ 
where z' = a 1 + iA' and A > A' > 0. 

The proof is analogous to that of Lemma 3.10 in [6] and will not be repeated here. 

We change (3 to 5 in the first denominator in the big bracket in (IA.22[) . the error is 
estimated by A 2 7?- 1 |a-/3||log7?| < 0{\ 1 / 2 ) using Lemma lA.21 Then we compute 



dQf(Q,y, 



a-Q + R- A 2 0(7) + ir] ql-Q + R- A 2 9(7) - it] 

f{Q,y,4>) 



i Im 



dQ 



a-Q + R- A 2 9(7) + irj ' 

Note that a + R > j + y 2 — X — 0(r) >A since g > 3A, so we can use the estimate 
(for e > e' > 0) 



Im 



x + ie 



-j dx = -tts(O) + Ole'/e) + Q{e'\ log e'| 



where the constants depend on H^Hc 1 - We obtain 
f(Q,y,4>) dQ 



Im 



o a - Q + R - A 2 G(7) + irj 



-Trf^a + R-X^i^^^^+Oiv/X) 

-vr/(7 + y 2 ,y,<A)+0(A) , 
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where 1Z is the real part of and in the last estimate we used the smoothness of / and 
\a — 7 1 = O(A). The result of these calculations is summarized as 

M 3 = — A 2 wr f j(7 + . 2 >3A)7(7 + . 2 ^^)d^ + y 

J a- f3 + 2G(j + y 2 ,0,y)-r-2i\ 2 l(-/)-2ir] 

as the contribution for Case 3. 

Using (|A,17j) and changing back the variables and adding the errors from Cases 1 
and 2, we obtain 

M = -Wit [ fAi + y 2 ,y,4>)J(i + y 2 ,y,ct>)dyd4> + 0(A i /2 ^ b A|) 

J a - [3 + 2G( 7 + y 2 , <j), y) ■ r - 2iA 2 J( 7 ) - 2irj 

__ox2- f f<rip)S(e{p) -7) dp , n / x l/2-9^ n 

" J a — /3 + 2(Ve)(p) • r - 2i[A 2 X( 7 ) + r?] + ^ A j ' 

A. 3 Proof of Lemma 17. 7L 

In this section we prefer to avoid the 2ir factors in the arguments of the trigonometric 
functions. Therefore we change variables p — ► we redefine the dispersion relation 

e(p):=e(|-) = ^(l-cos(p«)), (A. 23) 

the integration domain, T rf := [— vr,7r] d , and the triple norm, ||tt||P:= |||^p|||. The task 
of proving Lemma 17.71 with the redefined data is equivalent to the original formulation 
up to changing the universal constant. For simplicity, we then remove all tildes from 
the notation. Hence in this section the integration domain is p 6 [— it, ir] d and e(p) = 
£?=i(l - coe(pW)). As before, we work in d = 3 dimensions. 

To prove the estimates in Lemma 17.71 we first replace oj{p) with e(p) in the prop- 
agators at the expense of an extra factor A 2 r/ _1 = ^ _K /( 2 + K ) using a straightforward 
resolvent expansion: 

1 1 cA 2 1 + cAV 1 

<i ,,.., +; r~, m ; ; . : , < 



\a — u(jp) + irj\ \a — e(p) + irj\ \a — e(p) + irj\ \a — uj(p) + irj\ \a — e(p) + ir/\ 

(A.24) 

Therefore Lemma 17.71 will immediately follow from 



Lemma A. 3 



dp 1 cq 3 / 4 | lo. 



3 



I ■= / , , ^ - < ^ ,„' (A.25) 

|a — e(p) + zr/| |p — e(p + gj + ir]\ \lq\l 
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W ■= 1 1 ^ ■ r\ -nil 1 111 ^ Cr? " 7/ I'!, g??|3 (A-26) 

J \a - e(p) + ir]\ \p - e(p + q) + irj\ \jp - r\\\ |qj 

uniformly in r, a, (3, and 

^ da [ , ^ . n ^ , < C7? ' 1/ ,^° g??|3 . (A.27) 



id 



\a — e(p) + irj\ \a — e(p + q) + irj\ \q\ 



We remark that these estimates are far from being optimal, e.g. (|A.25|) and (|A.26|) 
are expected to hold with a prefactor r? -1 / 2 instead of r/ _3//4 and r?™ 7 / 8 . 

Before the proof of Lemma lA. 31 we establish the following 
Proposition A. 4 

/ 1 rfM -L ' I III 1 III ^ C ' ^^l 3 ' ( A ' 28 ) 

J \a — e{p) + ir]\ \jp — r\\\ 

! 1 7^ r-r -i 1 „, < crf K l 2 \ logr/l 3 . (A.29) 

J \a — u>(p) + ir]\ \jp — r\j 



Proof of Proposition A.J^, The second inequality follows immediately from the first 
one by (|A.24|) . For the proof of the first one, we consider the integral restricted on the 
set where — r||| = \p — r\ + rj, i.e. the closest critical point to p — r is 0. The other 
regimes are treated analogously. Since \p — r\ > \p\ — r\\, we have 

dp 1 f dpi f dp 2 dp 3 



\a — e(p) + irj\ \p — r\ + r\ J \p\ — r\\ + rj J \a — cospi — cosp2 — cosp3 + irj\ 
The inner integral is computed for any fixed p\ by using 



sup 

c 



f dp 2 dp 3 2 ( a ori\ 
— r < c log 77 . A.30) 

J \C — COSP2 — COS|?3 + IT]\ 



Then (|A.28|) follows after the dpi integration. 

To see (IA.30p . let f(p) := cosp2 + cosp3. On the set {(p2>P3) : \^f(p)\ < v)i we 
use the trivial bound on the integrand. The volume of this set is 0(rj 2 ), so this 
contribution is negligible. On the complement set, |V/| > rj, we use the coarea formula, 

1(|V/| > V)dp 2 dp 3 = f 2 da r 1(|V/| > y)d£ 
\c-f(p) + iv\ J-2 |c - a + irj\ / s(a) |V/(p)| 

with S(a) = {(p 2 ,P3) ■ f(p) = and dl is the arclength on S(a) that is an analytic 
curve with finite length. Since |V/(p)| > |||p||| + rj, the integral on the set S(a) is bounded 
by | log n\ since in the worst case it is a regularized l/\x\ singularity on a one dimensional 
piecewise smooth curve. □ 
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Proof of estimate iA . 26\) in Lemma \A.3[ We first show how the estimate ()A,26p 
follows from (|A.25|) . We distinguish two cases. If |||g||| < 100r ? 1 /8 ) then we estimate the 
(3 denominator trivially by rf~ 1 and use Proposition IA.41 The total estimate is 



ct] 1 \ logr/| 3 < 



ci] 7 / 8 1 log r/| 3 



Now we assume that I ? I > IOO77 1 / 8 . Let S := {p : \\\p - r|| < x^r/ 1 / 8 }. On the set 
S c we will estimate the |||p — J"!" 1 term by cry -1 / 8 and we use (|A.25p . 
Finally, we consider the integral on the set S. From 

inf (\ sinx\ + | sin(x + y)\) > — min{|y — kn\ : k G Z} , 

we see that 

|Ve(r)| + |Ve(r + g)| > ^|||g||| . 

Hence we can assume that |Ve(r)| > ^\\q\\ (the other case being identical). Since 
sup M ||e" (it) || < 1 for the usual matrix norm of the Hessian, we see that, for all p G S, 

|Ve(p)| > ^\\q\\ - r^V 1/8 > 4q\\ . (A.31) 

We again estimate the denominator trivially by r/~ 1 in J(r) and use the coarea 
formula to obtain 

! f l(peS)dp 1 _ _! f 6 da /• l(|b-r||| < T l 5 r ? 1 / 8 )da 



v l = V 1 

'a- e(p) + irj\ \\p - r\\ J \a - a + ir]\ J s{a) \\p - r\\ \Ve(p)\ 

where S a is the level surface {p : e(p) = a}, consisting of finitely many smooth pieces 
and da is the surface measure. Using (|A.3ip . we continue with 

! f l{p € S)dp 1 < or?' 1 /" 6 da /■ l(|b-r||| < T ^r ? 1 /8)da 



a - e(p) + ir?| |||p - r||| \\q\lJ l« ~ a + H J s(a) lb ~ r 

^ cr/~ 7 / 8 1 log 77 1 /" 6 do ^ cr/ -7 / 8 1 log 77 1 2 



/ |a — a + i7/| 1 9 1 

Here we used that the regularized l/|x| singularity on two dimensional regular surfaces 
E can be estimated as 

1(M < e)da(x) f l(\x\ < e 



This completes the proof of (|A.26p . □ 



f l(\x\ <e) ^ 
7r2 \x\+V 
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Proof of estimate iA . 25\) in Lemma \A.3l We fix q throughout the proof and we 
assume that |||g||| = \q\ + rj; the proof for the neighborhood of other critical points are 
analogous. We can assume that | 1 > 5 Ik II by permuting the indices and we can also 
assume q 2 G [— 7r/2, 7r/2], by a shift. 

We use the coarea formula 

dp 1 f 6 f 6 dadb f d£ 



I : 



a - e(p) + irj\ \/3 - e(p + q) + ir]\ J J \a - a + irj\ \/3 - b + irj\ J s ^ b) \F(p)\ ' 

(A.32) 

where we define the vector 

/ sinp 2 sin(p 3 + q 3 ) - sinp 3 sin(p 2 + 92) \ 
F{p) : = Ve(p) A Ve(p + q) = sinp 3 sin(pi + qi) - sinpi sin(p 3 + 93) , (A. 33) 

\sinpi sin(p 2 + 92) - sinp 2 sin(pi + q%) J 

the set 

S(a, b) := {p G T 3 : e{p) = a, e(p + q) = b) 

and let d£ denote the arclength measure. Since e(p) is analytic, the set S(a, b) consists 
of piecewise analytic curves with the exception of finitely many a, b. The lengths of the 
curves are bounded. 

Let f(p2) be the first component of F(p) as a function of P2 with fixed p 3 , 

/(p 2 ) := sinp 2 sin(p 3 + q 3 ) - sinp 3 sin(p 2 + 92) , (A. 34) 

and compute 

/'(p 2 ) ■= cos P2 sin(p 3 + q 3 ) - s'mp 3 cos(p 2 + 92) ■ 
For q 2 G [— 7r/2,7r/2] a simple estimate shows that 

\f{P2)\ 2 + \f'{P2)\ 2 > (| sm( P3 + g 3 )| + | sinp 3 |) 2 sm 2 f . 

Therefore 

\f(P2)\ + |/'(P2)| > cM sin(p 3 + ^3)1 + 1 sinp 3 |) \q 2 \ (A.35) 

uniformly in p 2 . We also note that 

\f"{P2)\ = I sin(p 3 +93)1 + 1 sinp 3 | . 

We now need the following one dimensional lemma: 

Lemma A. 5 Let f G C 2 [— vr,7r] satisfying 

inf(|/(p)| + |/'(p)|)><5, ||riU<M. 
p 

Then fore < 1/2 

VoZjp : |/(p)| < ed} < 8emin{87rM<r\l} . 
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Proof. We can clearly assume that eM < 5/32, otherwise the statement is trivial. 
Suppose now that for some po we have |/(po)| < sS. Then \f'(po)\ > 5/2. For defi- 
niteness, assume that f'(po) > 5/2. Clearly f'(p) > 5/A for all \p — po\ < \5M~ Y since 
f'{p) > f'{Po) — \\f"\\oo\p — Po\- For any p with 8e < \p — po\ < \5M~ X we thus have 
\f(p)\ > e5 and this set is non-empty since eM < 5/32. 

Define the intervals I\ := [p G [— vr, n] : \p — po\ < and I2 '■= {p S [ — vr, 7r] : 

\p — Pol < 8^}; then clearly I2 C I\ and \f(p)\ > s5 for any p E I\ \ Ii- Thus I2 occupies 
at most 

— < min{32Afe^ 1 ,4e/vr} 
\h\ 

proportion of I\. Repeating this argument for each interval where < 5e, we obtain 

the lemma. □ 

Using this lemma for / from (|A.34p . we obtain 

Vol{p 2 : |Fi(p)| < ce(\sm(p 3 + q 3 )\ + | sinp 3 |) |g 2 | } < ^ (A.36) 

for any e < \ and uniformly in all other variables. It is trivial to see that 

Vol|p 3 : I sin(p 3 + q 3 )\ + | sinp 3 | < p} < c/j 
for any \i > 0, uniformly in q 3 . Let 

Mj := {p 3 : 2~ j < I sin(p 3 + q 3 )\ + | sinp 3 | < 2- J+1 } 

then 

Vol(Mj) < c- 2- j . 
We consider these sets for j's such that 2~i > 2r7 3 / 4 and set 

M* := |p 3 : I sin(p 3 + q 3 )\ + | sinp 3 | < 2ri 3 / 4 } 

for the rest. 
Let 

M:= {( P2 ,P3) : |Fi(p)| < c^|g 2 |} 

then 

Vol(M) < Vol(M*) + Vol(M n M,-) 

3 

< cry 3 / 4 + ^ {(p 2 ,p 3 ) : P3 £ , |Fi(p)| < c£/i • 2 J '(| sin(p 3 + g 3 )| + | smp 3 \)\q 2 \} 

j 
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< c, 3 / 4 + V VoKM,)^ < c, 3 / 4 + C£/i , |lQg7?l 
~ ^ 5 92 ~ 92 

from (|A.36P with e modified to ep ■ 1? and using the uniformity of these estimates. At 
the end we will choose en. = r/ 3 / 4 , so the constraint 2~ J > 2r/ 3 / 4 will guarantee that the 
condition e < 1/2 in Lemma I A. 5 1 holds for the modified e. 
We split the I integral into two parts: 



I 



[1(M) + l(M c )]dp 1 
\a — e(p)+irj\ \f3 — e(p + q) + irj\ 



On M we use the trivial rj 1 bound for the second denominator and integrate out the 
other one in the p\ variable (independent of M) to bound its contribution to / by 

Cr ? -3/2( r? 3/4 + |logr?| .^ (A37) 

v |92K 



since 



c 



sup / ^1 < C7? -V2 . (A . 38) 

|c — cospi + vq\ 



On M c we use (|A.32p and the estimate \F\ 1 < (cefi\q2\) 1 to get 

c| log r/| 2 
e^l92| 



(A.39) 



The bounds in ()A.37p and (|A.39P are optimized for e/i = ry 3 / 4 and give (|A.25p , Finally 
the proof of (|A.27|> is almost identical, except that instead of (1A.38P we integrate out a 
to collect only | log r/| instead of ry^ 1 / 2 . □ 
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